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FOREWORD 


CHAM of North America Incorporated has performed a Rocket Injector Anomaly 
Study under the NASA Contract NAS3-23352. The purpose of the study was to modify, 
test and demonstrate a computer code for predicting three-dimensional two- 
phase spray flow and combustion in rocket engines. The modified computer code 
REFLAN 3D-SPRAY (REactive FLow ANalyzer 3-Dimensional, with two phase spray) 
and results of parametric studies have been described in the folloiwng two 
volumes: 

Volume 1: Description of the Mathematical Model and Solution Procedure; 

Volume 2: Results of Parametric Studies. 

Transfer of the code to NASA LeRC computer center, and preparation of a user's 
manual are recommended as next steps of the study. 

The authors wish to thank all those who have contributed to this work. In 
particular, thanks are due to Larry P. Cooper and Ken Davidian of the Communication 
and Propulsion Section of NASA LeRC; and to Laurence Keeton, Jack Keck, Kelli 
King, Janet Siersma, and Ronni Rossi c of CHAM NA. 
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SECTION 1 
SUMMARY 


The liquid fuelled rocket engine combustors consist of an injector plate and a 
thrust chamber. The injector plate consists of a number of propellant injectors 
which are designed to atomize the liquid jets of reactants and to promote 
intensive mixing between the vaporized components. 

Figure 1.1 shows the schematic of the rocket engine and injector plate, with 
L0X-RP1-L0X unlike triplet injectors, considered. The purpose of the study 
was to demonstrate an analytical capability to predict the effects of reactant 
injection non-uniformities (injection anomalies) on heat transfer to the walls. 

For this purpose an existing three-dimensional single-phase flow and combustion 
computer code (REFLAN3D: - REactive FLow ANalyzer, 3-Dimensional) has been modified 
for simulating two-phase flows in liquid propellent rocket engines. The modified 
code is referred to as REFLAN3D-SPRAY. 

Mathematical basis and solution procedure of REFLAN3D-SPRAY are described in 
Volume 1. This volume presents computational results to demonstrate the capability 
of the code. Reported results include numerical tests as well as several 
parametric test cases of the model rocket combustor (Figure 1.1). 

The results of the effort undertaken in this program can be summarized as 
follows: 

The REFLAN3D computer code with Eul erian-Lagrangean technique has 
been adapted for three-dimensional, elliptic, two-phase flow with 
evaporation, heat transfer and combustion in liquid fuel rocket 
engines. 
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- The code check-out calculations satisfied uniformity and symmetry 
requirements for both nonreactive and reactive turbulent flows. 

- The parametric studies of the two-phase flow and combustion in 
rocket engine combustor showed qualitatively correct response for 
variations of geometrical and physical parameters. 

- The injection nonuniformity test case, i.e. with blocked 25 % fuel 
injector holes near the chamber, showed significant changes in the 
central flame core and minor influence on the wall heat transfer 
(Figure 1 . 2 ). 

- The areas for model improvements, in both physical and computational 
aspects, have been identified and recommended for further studies. 


The recommendations for further studies include: 

a) preparation of a user's manual and code transfer to NASA LeRC; 

b) improvements in the physical models of evaporation, turbulent 
diffusion of droplets, and chemical reaction; 


c) numerical improvements to increase the accuracy of solution method; and 

d) verification studies with and without comparisons with experimental 
data. 


RO INJ ANAL* BASIC 21-XGRID 


XV PLANE 4 
TEMP CONTOURS 
PHIHIN 1 .S60E+02 
PH I MAX a . 525E+03 


CONTOUR LEUELS 

1 4.as0E+ea 

2 6.901E+02 

3 9.522E+02 

4 1.214E+03 

5 1.476E+03 
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7 2 . 000E+03 

8 2.263E+03 



a) full injection 
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b) blocked 25% of the fuel injector holes 



Figure 1,2 Temperature contours within the rocket, enginp 
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SECTION 2 
INTRODUCTION 

The objective of the present study was to demonstrate an analytical capability 
to predict the effects of reactant injection-nonuniformity (injection anomalies) 
upon local and overall heat transfer in a liquid propelled rocket engine 
combustion chamber. For this purpose the REFLAN3D (REactive FLow ANalyzer 
3-D imensionall ) computer code has been modified for the two-phase spray 
calculations arid applied in this study. The modified code is referred to 
as REFLAN3D-SPRAY. 

Detailed description of the mathematical model is provided in Volume 1, 
entitled "Description of the Mathematical Model arid Solution Procedure". 

This report constitutes the second volume of the documentation and describes 
the computational results and analysis of: 

a) code check-out test cases (section 3), and 

b) five parametric test cases (section 4). 

Based on the analysis of these results further model refinements are identified 
and described in section 5. The need for model varification against experimental 
data is also discussed in section 5. 
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SECTION 3 

CODE CHECK-OUT CALCULATIONS 


In order to predict the effects of injection nonuniformity on the fluid flow 
and heat transfer pattern the numerical model has to be capable of responding 
to small changes in injection pattern. Therefore, in addition to common checks 
of computational stability, convergence rate etc, a model accuracy for 
handling cyclic boundary conditions must be checked before it can be used 
for injection anomaly studies. For this purpose two test cases were selected 
to check solution symmetry and uniformity characteristics. 

3.1 SYMM E TRY TESTS 

Figure 3.1 presents the geometry configuration and computational grid used for 
the symmetry test calculation. The fluid enters the cylindrical duct through a 
partially blocked front plane. Two jets: 

- the inner jet with smaller velocity, and 

- the outer jet with larger velocity, 

enter through a V-shaped opening in the front plane. The fluid leaves the 
duct through fully open exit plane. 

A 180° sector of the duct has been used as the solution domain for the calculations 
with uniformly distributed grid (NX: NY: NZ - 21: 10:. 4). Due to the symmetry of 
the combustor geometry and symmetrical boundary conditions, the results of the 
computations should be symmetric about the 0=90° plane. 

Results of the following two test cases are presented in this section: 

A) A constant temperature (300°K), constant viscosity (.001 kg/ms) 

flow with fixed inlet velocities (10 and 20 -~) and fixed exit pressure 
(relative pressure, P .. = 0), specified as a boundary condition. 

tJX I l 

B) All conditions same as is case A, except that the larger inlet 
velocity is changed from 20 to 50 m/s. 

The k-e turbulence model has been used for both test cases. Figure 3.2 presents 
calculated velocity vectors in all four axial planes for test case A, Note 
that the minimum (VELMIN) and maximum (VELMAX) velocity vectors at symmetric 
planes (1 and 4) and (2 and 3) are identical. The shape and size of the 
recirculation zones on plane 1 and 4 are also the same. 
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XV PLANE 1 
UELOCITY PLOTS 
UELNIN 3« 160E—02 
UELNAX 9 . 825E+00 



XV PLANE 3 
UELOCITY PLOTS 


UELNIN 8 . 245E-03 


UELNAX 1.908E+01 



XV PLANE 3 
UELOCITY PLOTS 
UELNIN 8.E45E-03 
UELNAX 1 . 90 8E -+-01 



XY PLANE A 
UELOCITY PLOTS 
UELNIN 3 * 16 IE— 02 

UELNAX 9 ♦ 825E+00 



Fiqure 3.2 Velocity vectors for the symmetry test case. 
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To enable a thorough examination, results of test case B are presented in 

the form of the following tables: 

Table 3.1 - u-axial velocity (m/s) for each 0^ plane 

Table 3.2 - v-radial velocity (m/s) 

2 

Table 3.3 - rw-angular momentum (m /s) 

2 2 

Table 3.4 - k-kinetic energy of turbulence (m/s ) 

Table 3.5 - global convective fluxes through the chamber cross section (kg/s) 

Tables 3.1 to 3.4 indicate perfect symmetry about 0=90°. Results are 

symmetric at 0^ and 0^, and 0 ^ and 0^. Another symmetry verification is almost zero 

(10 ^) angular velocity at k=3 plane (0=90°). Note that "backward boomerang" 

for velocity staggering is employed. Further guidance on interpretation 

of the computer output is provided in Appendix A of the present report. 

Calculations were performed on a CRAY1 computer and total execution time for 
100 iterations was 5.6 seconds. Figure 3.3 presents the convergence rate for 
test case B in the form of residual error variation with number of iterations. 

The residual error e is calculated as the total mass inbalance for the entire 
calculation domain i.e. 

e = E ( EC d ) 
i j k d 

where i.ik denote grid indices, and d = N, S, E, W, H, and L signifies grid cells 
boundary at which the convective flux is calculated. 

The numerical convergence can also be verified from table 2.5 where global or 
total convective fluxes across various sections are presented. The notations 
used are: 

Global net convective flux 


FX. = 2 E (puA) k=l . . . . N=4; j=l . . . . M=10 

1 k j 

Global positive flux 

FXP. = E E max (0.,puA) 
k j 

Global negative (recirculating) flux 

FXM . = E E min (0. ,puA) 

k J 

where u - axial velocity, p - density and, A - grid cell face area, and 
i - axial, j - radial, k - circumferential coordinate indices. 
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Table 3.1 Axial Velocity 
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3.136 + 0 1 2,236 + 01 1 
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a, 556 + 01 2.08E + 01 1 
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3 ■■ -.Jiy, . , , ■ > 

-1,036+00 7,376-01 1 

3.63E+00 3.886+00 3 

8,a5E+00 7.Q26+00 5 

1 « 4 1 E + 0 1 1,096 + 01 8 

2,086+01 1.53E+01 1 

2,746+01 l,95g+01 1 

3,136+01 2,236+01 1 

" 37076 :+ 01 — 27266+01 1 

2,556+01 2,086+01 1 

2,186+01 1,916+01 1 


J Is t 
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r 9 O.OQE 
8 >< 0,09E 
7 0,006 
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'/•// 8 /;//■> ,11111 ’ 7 1/ 

•i,«3E+00 .3,486+00 
>2,436+00 “ 4,166*00 
3,476+00 4.826+00 

4,766+00 5.64E+0Q 

6,26E+0Q 6.57E+00 

7.88E+00 7.55E+00 

9.49E+00 8,«8E+00 

1 , 09E + 0 1 9.27E+00 

1 , 19E + 01 9.79E+00 

1.24E+01 1,00E+01 


(.47E+00 
1.98E+00 
> , 42E+00 
>,936+00 
>,50E+00 
\ 066 + 00 
r , 57E+00 
r ,97E + 00 
i,22E+00 
t , 30E + 0 0 


3, 18E+00 
3.88E+00 
4.45E+00 
5, 07E + 00 
5,736+00 
6,39E*0Q 
7 , 006 + 0Q 
7.53E+00 
7.95E+00 
8 , 22 E + 00 _ 

8 / > 

3, 18E + 00 
3,86E+O0 
4.456+00 
5,076+00 
5.73E+00 
6 i 396+00 
7 ,006+00 
7,~536+QO~ 
7.95E+00 
8,226+00 


'4,476+00 ? 4 
4 ',9 66+ 00-/5 
5,426+00 5 

5,9 JE+OQ 5 
6 .50E + 00 6 

7.06E+00 6 

7.57E+00 6 

7 , 97E + 00 . ,7 
8.22E+00 7 

0.3OE+OO 7 


», 916+00 
>,326+00 
> , 62E+00 
i, 956 + 00 
>,266+00 
>,60E+00 
>, 86E + 00 
’ , 05E + 0 0 
r, 15E + 00 
' , 176 + 00 


3.81E+00 
4.32E+00 
4.70E+00 
5,096*00 
5.50E+00 
5 , 9 1 E+00 
6,306+00 
6.65E+00 
5,936+00 
7,126+00 

. 9. ■ 

3.81E+00 
4, 32E + 00 
4 , 70E + 0 0 
5 , 09E + 0 0 
5.50E+00 

s o i 6+OQ 
5,306+00 
5765E+00 
6.93E+00 
7,126 +00 

; >■ -9 - 

4,9iE*oo: 

5+32E+0C 

5 , 62E+00 

5.95E+00 

b, 28E + 00 

5,606+00 

6 , 866+00 

7.05E+00 

7,156+00 

7,|7EtOO 


5,251+00 
5, 55E+00 
5,726+00 
5,886+00 
6,036+00 
6.16E+00 
6.25E+00 
6.30E+00 
6,326+00 
6.31E+00 


/_ U ) iLLLLLLL . 

, 35 E + 00 
, 7 1 E + 00 
, 94 E + 00 
, 15 E + 00 
. 37 E +00 " 

, 596+00 
, 8 1 E + 00 
, 0 1 E + 00 
, 17 E + 00 
, 28 E +00 > 

tt : j 

, 35 E + 00 
, 7 1 E + Q 0 
, 94 E +00 
i 15 E + 00 — 
. 37 E +00 
, 596 + 00 
, 81 E + 00 
1 0 1 E + 00 
, 17 E + 00 
, 286+00 i 


, 256+00 

, 556+00 

, 726+00 

.88E+00. 

i 03 E+ 00 * 

, 16 E + 00 

, 256+00 

. 30 E+ 00 ■ 

, 326+00 


CA/T 






**** CHECK2 PRINT OF UEB IlF 1 ItslO Jl* I JLs 4 

X3' 1 2 3 4 5 6 7 8 9 10 , 

4 5,776+00 5,926+Oe 5,926+00 5,856+00 5.78E+00 5.71E+00 I.67E+00 5.70E+00 5.84E+00 6,27E+00\ 

X= 3^5.326+00 5,066+00 4,066+00 4,946+00 4.93E+00 4.92E+00 4.93E+00 4.96E+00 5.02E+00 5,19E+00i 

21/5, 326+00 S, 066+00 8,966+00 4,946+00 4,936+00 4.92E+00 4.93E+00 4.96E+00 5.02E+00 5,196+00 

1 5,776+00 5,926+0? 5,926+00 5,856+00 5,786+00 5,716+00 5,676+00 5.70E+00 5,846+00 6,276+00^ 
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Table 3.2 Radial Velocity 


• v * **• it m $ w i v #’w> * it 9 ft ««*•*»*•* 9 


V-VEtOCITY ' 

* ** 


K <b 1 **ft 


? 9 

U i 


^ y.-r ' • •*>/' ? S* SSS' s r'r / , 

U 

c 

6 ,. y 


..... 8 • ■//,/* 

" • • // / >> t\ 

ioi 

-l,62£+0o 

-6.41E-01 

3.00E-01 

7. 29^-01 

' 4,792-61’" 

6, 48E-0 1 

“Tr782E»0l“ 

T TTH^r 

3.S1E-01 

4 , 60E-0 1 


--3.06E+0Q 

-1 , 02E + 00 

6, 22E-0 1 

l.aoE+oo 

1 , 59E+00 

1 , 23E+00 

9.312-01 

7,382-01 

6,842-01 

7.87E-01 

A 

-4.25E+00 

-1.20E+00 

9,892-01 

2.01E+00 

2 , 20E + 0 0 

1.77E+00 

1.35E+00 

1,082+00 

9.98E-01 

1,062+00 

A 

"5,OOE+Oo 

-1 ,15E + 00 

1 ,«12+00 

2.55E+0Q 

2,702+00 

2.22E+00 

1 ,722+00 

1.39E+00 

1 , 29E+00 

1.30E+00 

b'A 

-5.22E+00 

-8.49E-01 

1.862+00 

2.99E*00 

3.10E+00 

2,592+00 

2 , 03E + 00 

1,662+00 

1 # 54E+00 

} -52E + 00 

A 

-5.08E+00 

-2.58E-01 

2.31E+00 

3.31E+00 

3.37E+0O 

2.86E+00 

2,272+00 

1.882+00 

1.76E+00 

1,7 1E + 00 


-4,52E+Oo 

5.B8E-01 

2 , b82 + 00 

3,492+00 

3.50E+00 

2., 992 + 0 0 

2.41E+00 

2.03E+00 

1,912+00 

1 , 88E+00 

A 

-3,54E+Oo 

1 ,20E + 00 

2.64E+00 

3.43E+00 

3.40E+00 

2,932+00 

2.39E+00 

2.05E+00 

1 , 96E+00 

1.96E+00 

2 ’ 

-1.77E+00 

1.35E+00 

2,462+00 

2.85E+00 

2,802+00 

2, 4 1E+00 

1.99E+00 

1.74E+00 

1 , 68E+00 

1.73E+0O 

i) 

■ 0 | (H)S 

-9 t 002+-O4- 

-0-. (M)4hW)0 - 

■M92 + M- 

— 6-.O0E-+4O 

0.M6 + 40 

4U002 + 00 

4t«E + O0 

^-,002 + 00 

-OrOOE + OO- 


J 15*4 


?/,//* \ ///■■ * '// ' - /// 2 /// / / 1 

-i.Obt+Oo 1.05E+00 


10£ 

9^-l,48E+Oo 
e£-9,2UE-0l 
75| 1.17E+00 


2.S3E+00 
2,962+00 
5, 51E+00 
3,062+00 
1 .81E + 00 


r ..bA 5.asE+oo 

5 2.65E+00 

4 3,462=01 2,952-01 

3 -2 , 7 1 E+Oo -1.31E+00 
iNltl 2 -1.19E+0Q "1 ,51E+00 
1 -0...002+.00 


3 •< ■■■ 

“T,' 26 e-oT 

1 ,2bE+0G 
1 .bOE + OO 
1 ,b4E + 0C 
1 , 28E+00 
a , 3oE »0 1 
-6,99E-C1 
-1 ,7bE + 0(j 
=1,992+00 
— -e-r Q OE*^- 


• • -4 ■ • 5- - . . 

2.42E-01 1,842-02 


7 > 


4.21E-01 -1.25E-02 -1 
4.90E-CU -1.08E-01 “3 
3.76E-01 -3.08E-01 -5 
1.99E-02 -6.U9E-01 -6 
-b . 13E-0 1 -1,142+00 -1 
*1 ,42E+0C -1.71E+00 -1 
■2.13E+00 -2.19E+00 =2 
•2,252+00 -2,202+00 -1 
• UrMi^OO &t-£OE + CM3 Q 


■ , 1 , l u 


LLU. 


S 


i met 


■1 , O 6 E+O 0 1 , 0 3cVo 0 
9$Wl,48E+Oo 2.13E+00 

6y-9,2«E-0i 2 , 96E + 0 0 

7 A 1.17E+00 3.31E+00 

5.45E+00 3.Q6E+00 

5 2,652+00 1,612+00 

a 3.«6E-0i 2,95E-01 

3 -2.71E+00 -1.31E+00 
2 -1, 192+00 "1,S1E+00 
1 ■■-6-i- O - O E-aC^e 0 0E-+8-0- 


•3,- j 


i , 28E -o i 
1 ,262+00 
1 » 60E + 00 
l ,64E+oo 
1 .28E+00 
a , 362 =0 S 
■6,992-01 
•1.76E+00 
•1 .99E+00 


/< ■ 


2. 422-0 1' 
U.21E-01 
a , 90E-0 1 
3,762-01 
1.99E-02 
= 6 , 1 3 E -01 
■1 ,a2E + 00 
•2, 13E + 00 
■2.25E+00 


-0,00 £-+-&«• — U-, OOE + OO- 


1,842-02 -7 
-1.25E-02 -1 
-1.08E-01 -3 
-3,062-01 -5 
-6.49E-01 -8 
= 1,1 4E + 00 -1 
-1,712+00 -1 
-2,192+00 -2 
-2.20E+00 -1 
- O - rOOfr + S - O . 0 


J 7 ■»//!.> \ n > >r • 
1 Of- 1 ,0b 2 E + 0 0 
9 j -3 a 062+00 
8<-4.25£+Oo 
7 !-5,OOE + Oo 
bj-5. 222 + 00 
j?-5, OBE+Oo 
l>-4, 526+00 
17-3.54E+00 
V,- 1.77E+00 
i ff -0. OOE + O O 


/// - 2 * . 

-b.ajE-Ol 

•1 i 02E+00 

•1 ,20E+0Q 
"1,1 5E + 00 
-8.49E-01 
-2.58E-01 
5.2BE-01 
1.20E+00 
1,352+00 
0 ,< H )€ + 0-0 ■ 


,-3, 


O-LJLi-L' 


2.002-01 
b,22E"01 
9.89E-01 
1 .41E+00 
1 ,862+00 
2.31E+00 
2,682+00 
2,8aE+00 
2.46E+00 
0 , 0 - 0 2 + 0 0 


uJLllluu 


V.m-oi 

i ,aoE+oo 
2.01E+00 
2.55E+00 
2,992+00 
3.31E+00 
3.a9E+00 
3,a3E+00 
2.88E+00 
■8-j-OOE+M- 


b,n&i 


1 .59E + 00 
2 , 20E+00 
2.70E+00 
3. 10E + 00 
3.37E+00 
3,50E+00 
3,a0E+00 
2.80E+00 


,aaE-02 
.90E-01 
.52E-01 
.79E-01 
.84E-01 
,2eE+00 
, b7E + 00 
, 0 0E + 0 0 
, 9aE+00 
, OOE^-W) 

*** K 

,aaE-02 
, 9 0E«0 1 
.52E-01 
.79E-01 
,8aE»0 1 
, 2e>E + 0 0 
.67E+00 
. , 0 0E + 0 0 
.94E+00 
^O-OtVOO 
*** K 

JjLU 


1.26E-01 
2 , 7bE«0 1 ' 

а, 5a£«0 l 

б. 70E-01 
9.30E-01 
i e 22E+00 
1.53E+00 
1.76E+00 
l.bbEtOO 
■WJIE + Od- 

3 *** 

7 ^ 


Jl. 


_UL 


- U . 


■t.blE-01 -1, 
>3,3aE»01 «3, 
■5.23E-01 -5, 
■7.35E-01 -7, 
•9.70E-01 -1, 

■ 1 , 22E + 0 0 -1, 

■ 1 , U7E + 00 -1, 
■1.65E+00 -1, 
■1 ,52E + 00 »i , 
0,00 £.*.<10 — . 0 , 

• 8 " 


73E-01 
57E-01 
5SE-01 
72E-01 
01E+00 
25E+00 
48E + 00 
65E+00 
52£»00 
00£. * JQ 


7a8E-oi 
, 23E + 00 
.77E + 00 
.22E+00 
.59E + 00 
, 8bE + 0 0 
.99E+00 
.93E + 00 

.aiE+oo 

i<M)E *0 0 


1.26E-01 
-2.76E-01 
-4,5at-01 
«6,70£«01 
-9.30E-01 
-1.22E+00 
-1 .53E+00 
-1 ,7b£+00 
•1 .66E+00 
- OtOOE + 04 - 
s U *** 

■ /:' n / L / J , 1 


'1 , b 1E«0 1 »1 , 
•3.34E-0J -3, 
•5.23E-01 -5, 
• 7 , 35E»0l -7, 
■9.70E-01 -1, 
■ 1 . 22E + 00 • 1 , 
•1.47E+00 -1, 
-1.65E+00 -1, 
■i,52Ef00 -1, 
-O^OOE-MWl- 0, 


4.82E-01 
9.31E-01 
1.35E+00 
1.72E+00 
2.03E+00 
2, 27E+00 
2.41E+00 
2.39E+00 
1 , 99£* 0 0 


A! U-L 


8 


732-01 

57E-01 

5SE-01 

722-01 

012+00 

25E+00 

48E+00 

65E+00 

52E+0O 

OO E+ OO 

■•9 / '• t j 


3.78E-01 
7,382-01 
1,082+00 
1,392+00 
1 ,b6E + 0O 
1 , 88 E +00 
2 , 0 3E + 0 0 
2.0SE+00 
1.74E+00 


3. 

6. 

9, 

1. 

1. 

1. 

1. 

1. 

1. 


-fryO-OE + OO- 6-r0 0 E +-44 — 0, 


51E-01 
842-01 
982-01 
292+00 
542 + 00 
762+00 
912 + 00 
962 + 00 
682 + 00 
04 <+ »0 


-8 , 23E-02 
-2,402-01 
-4.36E-01 
-6.63E-01 
-9.15E-01 
-1,192+00 
•1 .46E+00 
-1 .67E + 00 
-1 .59E+00 
o^wifr+oo 


€*ir 

-k- 


-8.23E-02 
-2 , 40E-0 1 
-4.36E-01 
-6.63E-01 
-9,152-01 
-1, 19E+00 
*1 .46E+00 
-1 ,67E+00 
• 1 , 592+00 
- O - r O 0 E *- 0 U ) — 

■ • / lo---.' 


€»r 

--$r 



4.60E-Q1 
7.B7E-01 
1.06E+00 
1.30E+00 
1.52E+00 
1.71E+00 
1.88E+00 
1 , 96E+00 
1.73E+00 
-& T - ME +- 04 — - 


t&r 


**** CHICK2 fRINT OF VNB II* i IL*15 Jl* i JL* 4 “ 

J I* « 5 4 5 

a 0,00f*0e »,0*E+0» . 0, OOE+OO 0, 00E+00 0, OOE+OO 

T“0,O9B+0r' »,0»E+0»“ ^,00E+00 0.00E+00 0, OOE+OO 

2 0,001+06 9,801+09 0,002*00 0, OOE+OO 0, OOE+OO 

1 0.00g+08 i.OOE+OS 0, OOE+OO 0,092+00 0, OOE+OO 


8 7 • 9 10 

0, OOE+OO 0. OOE+OO 0, OOE+OO 0, OOE+OO . 0, OOE+OO 

0, OOE+OO 0^ OOE+OO 0 , OOE+OO ' 0, OOE+OO' O.OOE+OO " 

0, OOE + OO 0, OOE+OO 0, OOE+OO O,OOE+O0 0, OOE+OO 

OgOOE+OO 0, OOE+OO 0,901+00 0, OOE+OO 0.00E+00 



Table 3.3 


Angular Momentum 
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V. 
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0 , 0 ot ♦ 0 0 
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O.OOE+Oo 

0,006+00 

O.OQE+OO 

O.OOE+OO 

O.OOE+OO 

0, OOE + OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

0, OOE + OO 

8 

7 

/ 

O.OOE+Oo 

O.OqE+09 

O.OOE+OO 

O.OOE+OO 

0.00E+00 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

1 

2 

0,006+0o 

O.OOE+OO 

O.OOE+OO 

0 v OCE + OO 

0 • COt + GO 

G , GutrCG 

G , G CE*0 G 

GtCOE^OO 

0| wQt^vO 

0 , OCE + OO 

8 3 

0 , OOE+Oq 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

0 , OOE + OO 

0, OOE+OO 

0, OOE + OO 

0. OOE + OO 

5 

A 

0,006+00 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

■j- 


O.OOE+Oo 

0 , OOE + OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 



O.OOE+OO 

0.GSE+00 

C 9 OOE+CC 

0 1 0 OE + OG 

O.OOE+OO 

0 , OOE+OO 

0,0OE+00 

C a OOE + OO 

O.OOE+OO 

O.OOE+OO 


A 

0,006+00 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

o.ooE+Oo 

O.OOE+OO 


1 

O.OOE+Oo 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

0. OOE + OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 
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T 3 ,//// 

>/,■: 

' ' • //3// /////." / 4 ,, / , 

•• • , 5 •//,. 

• > j 6 //_*/; 

> >///' 7 > .•! 8 '///. 

• '/ / 9 /t .■ 

•'•".10 X " 


's 

1 

•e! • 8$E + 00 

-5. 02E+OO 

-i .296+00 

6 , 4 6 E - 0 1 

1 ,916 + 00 

2,306+00 

2.36E+00 

2.32E+00 

2.07E+00 


i 

2.19E+00 

-1 , 76E+00 

-2.57E+00 

-1 .55E + 00 

-7.21E-02 

l.OOE+OO 

1.45E+00 

1 .5BE + 00 

1.58E+00 

1.37E+00 

c 


2.86E+00 

-8.34E-01 

•2 • 07E+00 

-1.55E+00 

-4 . 93E-0 1 

3.62E-01 

8.13E-01 

9 1 69E«0 l 

9,91 £-0 1 

8.31E-01 

i 


3,28fc + 0[) 

-i .eiE-oi 

-1.62E+00 

-i ,466 + 00 

-7.53E-01 

-6,966-02 

3, 19E-01 

4.80E-01 

5.12E-01 

4.02E-01 


✓I 

3.11E+00 

1.67E-01 

-1 .24E + 00 

-i. 326+00 

-8.77E-01 

-3.71E-01 

-4.46E-02 

1.04E-01 

1.39E-01 

6.82E-02 

Z; 


2.01E+00 

2.T3E-01 

-9.28E-01 

-1..13E + 00 

-8.92E-01 

-5 , 42E-0 1 

-2.87E-01 

-1 , 62E-0 1 

-1.31E-01 

-1.73E-01 



1 .01E + 00 

1.20E-01 

-7.23E-01 

-9,366-01 

-8.24E-01 

-6.00E-01 

-4, leE-01 

-3.20E-01 

-2.95E-01 

-3.20E-01 

, v 


3,l?E-0j 

-1.37E-01 

•5 . B8E-0 i 

-7, 32E-01 

-6.87E-01 

-5.57E-01 

-4.35E-01 

*3 a 68E»0 1 

-3.51E-01 

-3 . 66E-0 t 

c 


lalOE-Oj 

-2.35E-01 

-4. 15E-01 

-4.65E-01 

-4.72E-01 

-4.07E-01 

*3 , 37E-0 1 

-2.97E-01 

-2.87E-01 

•2 , 98E-0 1 


- 

5 , 186 = 03 

-6,976=02 

-1 .08E-01 

- 1 , 23E-0 1 

-1 .21E-01 

-1,046-01 

-8.74E-02 

-7,8 iE-02 

-7,676-02 

-8 , 07E-02 
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} -8.12E.08 
/ 6.86E-08 
J.4.68E-08 
•^,3, ORE-Os 
i -1.73E-08 
,1.926-08 
1.686-08 
1 ,876-0? 
s-lie-OA 


”5 »2bE-0S 
-4.03E-06 
-S.43E-06 
-2.98E-06 
-2 . 5 1 E-0 6 
-2.iBE.08 
-2.17E-0* 
-1 .69E-08 
-1.61E-09 
1 .57E-0S 


■3.34E-06 
•2.58E-06 
-2.12E-0S 
■ 1 , 86E-08 
-1.67E-0B 
■1.52E-08 
»1 1 37E-08 
■9.59E-09 
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Table 3.4 Kinetic Energy of Turbulence 


K.E, OF TURBULENCE 
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0,0000006+00 

FULFLXa 

0 , OOOOOOE+O 

1 % 

8 

XMa 
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FXa 
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Table 3.5 Convective Fluxes Through the Chamber Crossections 
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The global mass flux FX = 14.73138 should be, and is preserved, at each 

axial location 1 = 1 1CL 

3 . 2 UNIFORMITY TESTS 

The uniformity tests have been performed on similar chamber geometry with axi- 
symmetric fluid entry (figure 3.4). A turbulent flow with propane-air combustion 
in a chamber with sudden enlargment has been considered. Figure 3.4 presents 
the geometry and inlet conditions. Fuel enters the chamber through the axial 
slot at x=0 with u=10^, v=0, w=0, and air through the annular inlet with 
u=50^, v=0, w=o (see table 3.6). Results are presented by the following tables: 

3.6 - u -axial velocity (m/s) 

3.7 - v -radial velocity (m/s) 

2 

3.8 - rw-angular velocity (m '/s) 

2 2 

3.9 - k -kinetic enerty of turbulence (m /s ) 

3.10- T-absolute temperature (°K) 

3.11- -fuel mass fraction (-) 

3.12- FX, FXP, FXN and FULFLX fluxes (kg/s). 

Calculated results are uniform on all 0^ (k = 1, 2, 3, 4,) planes. Note that 

the recirculation zone extends up to the middle of the chamber. The angular 

-8 -12 

momentum rw should be zero throughout the chamber. It was predicted as 10 - 10 

at 0=90° and as 1 0~ 4 - 10' 6 at 0 = 45° and 0 = 135° respectively. These values are 
within acceptable error limits. 

The values k, T and also indicate good uniformity. Table 3.12, in 
addition to global mass fluxes (FX, FXP, FXN) S also provides total cross- 
section unburned fuel flow rate, which is defined as: 

FULFLX = Z 5 (pu«.m fu ) J>k 

The convergence and conservativity of the numerical algorithm can be verified 
by FX distribution. Note - the FX : = 17.79037 kg/s value is exactly preserved 
at each axial i-plane (table 3.12). 

The convergence rate for the uniformity test case is presented in figure 3.3. 

The residual error for test case B indicates better convergence for the flow 
with uniform inflow. Intensive recirculation zones in case A (symmetry test) 
are responsible for relatively slower convergence rate. 


3-11 



3-12 



12 3456 789 10 

ENTRY PLANE 

Figure 3.4 Combustion Chamber Geometry 
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Table 3.6 Axial Velocity 
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3.77E+01 

3.816+01 

3,866+01 



5 s 006*0 J 

5 s 97E + 0 1 

3.29E+01 

2 , 75E + 0 1 
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4.61E+01 

4.276+01 

3.78E+01 

3.22E+01 

2.51E+01 

4 ,786 + 01 

4.29E+01 

4,206+01 

4,176+01 


riv"^ 

i , 00t*Q 1 

3, 18E+01 

3,286+01 

3 , 20E + 0 1 

2 , 96E + 0 1 

2,346+01 

4.81E+01 

4.35E+01 

4.26E+01 

4,236+01 


FU i 

1 , QQE+0 1 

2.17E+01 

2,576+01 

2.69E+01 

2.65E+01 

2.13E+01 

4,766+01 

4.34E+01 

4.27E+01 

4.26E+01 > 









*** K 

s 2 *** 





j 

I* 1 


2 

3 

4 

5 

6 

. 7 

8 

9 

10 


10 

0,006 

+ 00 

-8.49E+00 

• 4 , 60E + 00 

-4.97E+00 

-1.90E+00 

3.96E+00 

2.83E+01 

3.26E+01 

3.346+01 

3.37E+01 


9 

0,006 

+ 00 

"5 , 3 1E+00 

"2.86E+00 

4,216-02 

2.75E+00 

7 , 45E + 00 

2.92E + 0 1 

3.37E+01 

3, 49E + 0 1 

3.57E+01 


8 

O.OOE 

♦ oo 

-3.75E-01 

3, 13E+00 

5 , OOE + OO 

6.91E+00 

1.05E+01 

3,096+01 

3,466+01 

3.58E+01 

3,676+01 


7 

O.OOE 

+ 00 

8, l&E+OO 

1.08E+01 

1 ,126+01 

1.19E+01 

1 , 38E+01 

3.37E+Q 1 

3,596+01 

3,696+01 

3.76E+01 


6 

0.006 

+ 00 

2,136+01 

2 , OBE + O 1 

1 .S6E + 01 

1 ,766 + 01 

1.74E+01 

3,746+01 

3,776+01 

3,816+01 

3, 8*6+01 


5 

5,OOE*Oi 

3.97E+01 

3.29E+01 

2.75E+01 

2.39E+01 

2, 09E+0 1 

ft, 16E+01 

3,976+01 

3,956+01 

•3V97E+01 


AlK* 

5,006*01 

4.62E+01 

ft , 12E + 0 1 

3.50E+01 

2.95E+01 

2.38E+01 

«,556+01 - 

ft,15E+01 

ft, 096+01 

“ft, 086*01" 


M 3 

5,006 + 0 t 

4,616+01 

4.27E+01 

3.78E+01 

3.22E+01 

2.51E+01 

ft, 786 + 01 

a,29E+01 

4.20E+01 

4.17E+01 


_ T 

1 ,006+01 

3.136+01 

3,286+01 

3 , 20E + Q 1 

2 , 96E + 0 1 

2.34E+01 

ft,81E+01 

ft, 356 + 01 

ft , 26E + 0 1 

4,236*01 


vu i 

l.OOE+Ol 

2.17E+01 

2.57E+01 

2 , 69E+0 1 

2.65E+01 

2.13E+01 

ft, 766 + 01 

«,3ft6+0i 

4,276*01 

4,266*01 j. 









■kick }\ 

5 j 



* 


J 

Is 1 


3 

3 

a 

5 

6 

7 

6 

9 

10 


so 

0,006 

+ST 

"'*g,S96 + do ' 

-S', 60E + 00 

• 4 # 9 7 ^ Q~Q 

-r,9oe*QQ~ 

' 3 , 966 + 00 

2,836 + 01 - 

3,266+01 

3.34E+01 

3;37E*ai 


9 

O.OOE 

+ 00 

-5.31E+O0 

-2.36E+00 

4.21E-02 

2.75E+00 

7 , 45E + 00 

2.92E + 0 l 

3 , 37E + 0 1 

3.49E+01 

3.S7E+01 


8 

0,006 

♦ oo 

-3.75E-01 

3,136+00 

5.0QE+00 

6.91E+00 

1 , 05E+0 1 

3,096+01 

3,466+01 

3 , 58E + Q 1 

3.67E+01 


7 

O.OOE 

+ 00 

6.16E+00 

1.08E+0! 

1, 12E + 01 

1, 196 + 01 

1.38E+01 

3.37E+01 

3.596+01 

3,696+01 

3,766+01 


6 

O.OOE 

+ 00 

2.13E+01 

2,086+01 

1.88E+01 

1,766+01 

1 , 74E+0 1 

3.74E+01 

3,776+01 

3.81E+01 

3,866*01 


T 

5.00E+01 

3.97E+01 

3.29E+01 

2.75E+01 

2.39E+01 

2 , Q9E + 0 1 

4,166+01 

3.97E+01 

3,956+01 

3 , 97E + 0 1 



5.00E+01 

4 , 62E + Q 1 

4,126+01 

3.50E+01 

2.95E+01 

2,386+01 

4,556+01 

ft , 15E + 0 1 

4.09E+01 

4,086+01 


3 

5 . QOE + O 1 

4.61E+01 

4.27E+01 

3 . 78E + 0 1 

3.22E+01 

2,516+01 

4.78E+01 

4,296+01 

4,206+01 

4,176+01 


“? 

l.OOE+Ol 

3.1SE+01 

3 ■ 286 + 0 1 

3,206+01 

2.96E+01 

2.34E+01 

4,816+01 

ft, 356 + 01 

ft , 26E + 0 1 

4,236+Oi 


pu i 

l ,OOE + Ol 

2,176+01 

2.S7E+01 

' 2.69E+01 

2 , 65E + 0 1 

2, 1 3E+0 1 

4,766+01 

ft,3ftE+01 

4,276+01 

4,266+01 i 









*»* 7C 

3 Hi** 



■ - ■ — — <$r 


j 

Is 1 


2 

3 

4 

. 5 ... . 

6 

7 

8 

9 .. 

10 _ 


10 

0,006 

+ 00 

-8.U9E+00 

-8.59E+00 

•4.97E+00 

-1.89E+00 

3.96E+00 

2,836+01 

3,266+01 

3,3«E+Q1 

3.37E+0I 


9 

O.OOE 

+ 09 

"5 , 3 1 E + 0 6 

•2 # 86E+00 

4,576-02 

2.75E+00 

7 ,ft56+00 

2,926+01 

3.37E+01 

3 , 49E + 0 1 

3.57E+01 


8 

O.OOE 

+ 00 

-3,73E-0i 

3,136+00 

5,016+00 

6,916+00 

1,056+01 

3.096+01 

3 , ft6E + 0 1 

3.58E+01 

3.67E+01 


7 

0,006 

+ 00 

8 , 16E + 00 

1 , 08E + 0 1 

1.12E+01 

l , 19E + 0 1 

1,386+01 

3,376+01 

3.S9E+01 

3 , 69E + 0 1 

' 3.76E+01 


6 

O.OOE 

+ 00 

2.13E+01 

2.Q8E+01 

1,686+01 

1.76E+01 

1.74E+01 

3.74E+01 

3 , 77E + 0 1 

3.81E+01 

3.86E+01 


n 5 

5,006+01 

3,«7E+0t 

3.29E+01 

2.75E+01 

2.396+01 

2,096+01 

4.16E+01 

3.97E+01 

3,956+01 

3.97E+0 1 


A\(U 

5, OOE+O 1 

4,626+01 

ft, 12E+Q1 

3,506+01 

2.95E+01 

2.38E+01 

4,556+01 

ft , 15E + 0 1 

4,096+01 

4,086+01 



3 

5.001*01 

4,616+01 

4.27E+01 

3 ,786+01 

3.22E+01 

2,516+01 

4,786+Oi 

ft , 29E + 0 1 

4,206+01 

4.17E+01 



1,006+01 

.-3.16E+0J 

J.28E+01 

-,3,20E+0l 

-2.96E+01 

3.34E+01 

«,81E+01 

ft , 356+0 1 

ft, 266 + 01 

ft,23E+01 



T.OOE+01 

* J.lfS+OT 

- a.sn+or 

"2.69E+01 

2.65E+01 

2, 13E+01 

4,766+01 

«;3«E+01 

~ft,27f+0r 

~a;26E+oi ~r 




, ,Vr 

- -V n • 





■ ^ ."v 





. . 


■ “ 


* 










fSjiwT rif*- . r 

TLplB "JJ* 

■ i Jte ft 



. . ■ ■■->- 




j 

IS- * 

'a--* 

^' 3 

- 'ft 

■ ' s 

6 

7 

8 

9 

10 


4 

4,aa68E + 91 4.1981E + 01 — ft,1368E+Q 1 - ft. 

0525E + 0 1 3 

.9491E+01 

1,846*6+01 

3, 75946+01.3, 69266*01 3,64096*01 

3,530*E+«r:~ 

3 

ft, 22606+91 *,19686+61 ft, 1389E+01 ‘ ft. 

0527E+01 3 

;9«936+01 

I.8«70E+01 

, '3i7599E*01—3;69Jir+0l—rj, 64156+01“ 

-3,53136*06—" 

2 

4.2269E+61 4, 1968E+01 0.1389E+01 4, 

05276+01 3 

,949 JE+Oi 

3,8«70E+0! 

3,75996+01 3,69316+0 1 3,6«15£+01 

3, 53136+01 


1 *+ 2?»0e+0y M98^E*81- -0,1 36 86+01 ft, 0 525 £ tO 1 5,949 16+0 1 5 , 8«**€+0 1 3, 7594Et01 5,& 9 2tC»OJl 3.6S09E+01 j,5306E»0| 
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Table 3.7 Radial Velocity 


• *«**»**±^#w*»*«*»*»*»*»*w V^VELOCXTT •* 

*** K m 1 *** 


J I 

= 1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

to 

• 1 s ?4g*0(V 

-2 . isr-na 

7,445-01 

6,495-91 

1.18E+00 

4.Q3E+00 

6.72E-01 

5.31E-02 

-3.66E-02 

2.50E-0 1 

9 

-3, 046+00 

4.69E-01 

1 ,412 + 00 

1,262+00 

2.16E+00 

7 , 32E+00 

1.33E+00 

1.73E-01 

4 ,692-04 

3.19E-01 

a 

“3 , 51E + 0 o 

1.25E+00 

1.98E+00 

1 , 79E + Q0 

2.95E+00 

1.02E+01 

1 .90E + 00 

3,016-01 

5 , 66E-02 

3.14E-01 

7 

•2.36E+00 

2.01E+00 

2.38E+00 

2. 16E + 00 

3.45E+00 

1.24E+01 

2 ,322 + 00 

4.08E-01 

1.14E-01 

2.79E-01 

t> 

t .676+00 

2.29E+00 

2.43E+00 

2.26E+00 

3.57E+00 

1,356+01 

2.50E+00 

4 , 77E-0 1 

1.71E-01 

2.50E-01 

5 

”2, 73£ = 0 i 

1 .2SE + 00 

1 ,812 + 00 

1 .94E + 00 

3.23E+00 

1.24E+01 

2.2SE+00 

4.89E-01 

2.27E-01 

2.S1E-01 

4 

-1.24E+00 

4.87E-01 

9.23E-01 

1 .21E + 00 

2.41E+00 

8.S8E+00 

1 .80E + 00 

4.77E-01 

2,602-01 

2.43E-01 

3 

»3,10E+0 0 

-2.62E-01 

9.34E-02 

3.80E-01 

1.38E+00 

4.12E+00 

1 , 2BE+00 

3.97E-01 

2.30E-01 

2.00E-01 

2 

-1. OOE+OO 

-3,106-01 

-7,245-02 

6.82E-02 

6 , 06E-0 1 

l , 32E+00 

6 , 69E-0 1 

2.31E-01 

1.37E-01 

1.17E-01 

1 

0 , OOEtOo 

0, OOE + OO 

0, OOE+OO 

0, OOE+OO 

0, OOE+OO 

0, OOE+OO 

0, OOE + OO 

0, OOE + OO 

0, OOE + OO 

0,006+00 







*** K 

5 2 *** 




J I 

= 1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

to 

-1 .74E + 00 

-2 , 19E-02 

7,445-01 

-6,492-01 

1 , 18E+00 

4, 03E+00 

6,732-01 

5.34E-02 

• 3 , 63E-02 

2.50E-01 

9 

-3, 046 + 00 

4.69E-01 

1 .41E + 00 

1 , 26E + 00 

2, IbE+OO 

7 , 33E+00 

1.33E+00 

1.74E-Q1 

1.11E-03 

3.20E-01 

8 

•3.53E+0Q 

1 .25E+00 

1,982+00 

l , 79E+00 

2.95E+00 

1 , 02E + 0 1 

1 .90E + 00 

3.02E-01 

5,762-02 

3.15E-01 

7 

•2 , 36E+QQ 

2.01E+00 

2.38E+00 

2, 16E + 0Q 

3,456+00 

1.24E+01 

2 , 32E + 0 0 

4,096-01 

1 , 16E-01 

2.80E-01 

6 

1,676+00 

2.29E+00 

2.43E+00 

2.26E+00 

3,586+00 

1 .35E+01 

2.50E+00 

4,792-01 

1 .73E-01 

2.51E-01 

5 

-2,736-01 

1.28E+00 

1.81E+00 

1.94E+00 

3,236+00 

1.24E+01 

2.28E+00 

4 , 9 1E-0 1 

2.30E-0 1 

2.53E-01 

4 

-1.24E+00 

4 , 87E-0 1 

9,242-01 

1 ,21E + 00 

2.41E+00 

8.58E+00 

l .80E + 00 

4 , 80E-0 1 

2,622-0 1 

2.45E-01 

3 

-3, 10E+00 

-2.62E-01 

9,3 7E-02 

3.81E-01 

1,386+00 

4,122+00 

1 .28E + 00 

4 , OOE-O 1 

2.33E-01 

2.03E-01 

2 

-1,006+00 

-3.10E-01 

-7 , 22E-02 

6 , 86E-02 

6,076-01 

1 .33E+00 

6,72E«01 

2.35E-01 

1.40E-01 

1,196-01 

1 

0 , OOE+OO 

0, OOE + OO 

0, OOE + OO 

0, OOE + OO 

0 , QOE + OO 

0,00£+00 

0, OOE + OO 

0, OOE + OO 

Q, OOE + OO 

0, OOE+OO 







*** K 

s 3 *** 




J I 

= 1 

2 

3 

4 

5 

b 

7 

8 

9 

10 

10 

-1,746+00 

-2.19E-02 

7.44E-01 

6.49E-01 

1.18E+00 

4 , 0 3E + 0 0 

6.73E-01 

5.34E-02 

-3.63E-02 

2.50E-QI 

9 

-3, 04E+0Q 

4.69E-01 

1.41E+00 

1 .26E+00 

2, 16E + 00 

7 , 33E+0U 

1.33E+00 

1,746-01 

1.11E-03 

3.20E-01 

6 

-3 ,S3E + Oo 

1,252+00 

1.98E+00 

1 ,792+00 

2.95E+00 

1.02E+01 

1 .9OE + 00 

3.02E-01 

5,762-02 

3, 15E-01 

7 

*2 , 36E+00 

2,012+00 

2.38E+00 

2 , 162 + 0 0 

3,456+00 

1 ,242+01 

2.32E+00 

4 , 09E-0 1 

1.16E-01 

2.80E-01 

6 

1 ,67E+0a 

2,296+00 

2.43E+00 

2.26E+00 

3 , 58E + 00 

1 .35E+01 

2.50E+00 

4.79E-01 

1,736-01 

2.51E-01 

5 

-2 , 73E-0 i 

1 ,23E+00 

1 .81E+00 

1,942+00 

3,236+00 

1.24E+01 

2.28E+00 

4.91E-01 

2,306-01 

2.53E-01 

4 

-1 ,246+00 

4 ,8?£-0 t 

9 - 24E-0 1 

1,212+00 

2,416+00 

8.58E+00 

1 ,80E + 00 

4,802-01 

2 , 62E-Q 1 

2 , 45E-0 1 

3 

*3, 10E+00 

-2.62E-0 i 

9.37E-02 

3,812-01 

1 .38E + 00 

4, 12E+00 

1,286+00 

4 , OOE-O 1 

2.33E-01 

2 , 03E-0 1 

2 

-1, OOE+Oo 

•3,102-01 

-1 , 22E-0 2 

6.86E-02 

6.07E-01 

1.33E+00 

6,722-01 

356-01 

1.402-01 

1.19E-Q1 

1 

O.OOE+Oq 

9,006+00 

0, OOE + OO 

0, OOE + OO 

0, OOE + OO 

0, OOE+OO 

0, OOE + OO 

0, OOE + OO 

0 , OOE + OO 

0, OOE + OO 







»** K 

* 4 *** 




J I 

= 1 

a 

3 ~ 

4 

' ' : 5 

6 

7 

6 

9 

10 

to 

-1,746+00 

-a,18E-02 

7.44E-01 

6,492*01 

1 , 18E + 00 

4, 03E+00 

6.72E-01 

5.31E-02 

-J,66E-02. 

2.50E-01 

9 

-3,o«e+oo 

4, 696-01 

l ,412+00 

1.26E+Q0 

2.16E+00 

7.32E+00 

1.33E+00 

1.73E-01 

4.69E-04 

3.196*01 

8 

»3,5}E+00 

1,256+00 

1.90E+00 

'1.79E+00 

*‘2.956+00 

" 1 ,086+01 

1.90E+00 

r.o iE-oi 

' S.66E-02 

3;i46*0l 

7 

•2,366+00 

S.01E+00 

2.J85+00 

2 , 16E+00 

3.45E+00 

1.24E+QI 

2 , 32E+00 

4.08E-01 

1 , 14E-0 1 

2, 796*01 

6 

1,616+00 

S, 296+00 

2.43E+00 

2.26E+00 

3.57E+00 

1,35E+01 

a,50E+00 

4,776-01 

1, 716*01 

2 , 50E-0 1 

5 

-2.736-01 

1.26E+06 

1.816+00 

1.94E+00 

""3,236+00 

l ,24E+01 

2.28E+00 

4.89E-01 

2.27E-01 

2.51E-01 

4 

-1 ,24E+00 

4,876-01 

9.23E-01 

1.21E+Q0 

2.41E+00 

8.58E+00 

1.80E+00 

4.77E-01 

2 , 60E-0 1 

2.43E-01 

3 

-3,106+00 

-2.6JE-01 

9.34E-08 

3.80E-0 1 

1 ,38E + 00 

4.12E+00 

1 , 28E+00 

3.97E-01 

2,306-01 

2.00E-01 

2 

•1,006+00 

-3,106-01 

-T.«6+oa 

6.82E-02 

6,062-01 

1.32E+00 

’ 6,69E*01 

2.31E-01 

T, 376-0 1 * 

1.17E-01 

1 

0,006+00 

0, OOE + OO 

0, OOE + OO 

0, OOE+OO 

0, OOE + OO 

0, OOE+OO 

0, OOE + OO 

0, OOE+OO 

0, OOE + OO 

0, 006+00 

*** CH6CK2 ?SJNT Of VN6 • -I H i 

11,610 Jl» 

1 Jt« 4 







1 

a 

3 

a . 

5 

t 

7 

8 

< 9 

.10 


a O.OOOOE+OO d f OOO0EtOO O.OOOOE+OO "O.OOOOE+OQ O.OOOOE+OO G.000QE+0Q O.OOOOE+OO ~0,0000E+00 0, OOOOE+OQ O.OOOOt+OQ 


3 

0 , 00002+00 

O.QOOOE+OO 

0 , OOOQE+OQ 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0,00006+00 

O.OOOOE+OO 

0,00006+00 

o.ooooe+oo 

2 

O.OOOOE+OQ 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0,00006+00 

O.OOOOE+OO 

O.OOOOE+OO 

0.00006+00 

-1 

0 . 0 

D'O 

o 

6 

o>o 

0.0 

0.0 

0.0 

0.0 

0\O 

0,0 
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SWIRL 


J 1 

= 1 

2 

3 

4 

5 

10 

0 ,006+00 

0, OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

9 

0 , OOE+Oo 

o.ooe+oo 

O.OOE+OO 

0.006+00 

0, OOE + OO 

8 

0 , OOE + Oo 

0 , OOE + OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

7 

0 , OQE+Oo 

0,006+00 

O.OOE+OO 

0 ,006 + 00 

O.OOE+OO 

6 

0, OOE + Oo 

O.OOE+OO 

O.OOE+OO 

0 .OOE + OO 

0, OOE + OO 

5 

0 , OOE+OO 

0, OQE + OO 

0 ,006+00 

O.OOE+OO 

0,006+00 

4 

0, OOE + Oo 

O.OOE+OO 

0,006+00 

0,006+00 

O.OOE+OO 

3 

0, OOE+Oo 

O.OOE+OQ 

0 , OOE + OO 

0,006+00 

O.OOE+OO 

2 

O.QQE+Oo 

0, OOE + OO 

O.OOE+OO 

0,006+00 

0,006+00 

1 

0, OOE+Oo 

o.ooe+oo 

O.OOE+OO 

0 .OOE+OO 

O.OOE+OO 

l I 

= 1 

2 

3 

4 

5 

.0 

-6,776-04 

-2.45E-04 

-7 , 18E-04 

-6 , 62E-04 

-1,776-04 

9 

-1,086-03 

-2.74E-04 

-5,946-04 

-5 , 27E-04 

-1,026-06 

8 

-1.276-03 

-3.23E-04 

-4,486-04 

-3 , 84E-04 

8.00E-05 

7 

-1,246-03 

-3.99E-04 

-3.21E-04 

-2,496-04 

1 ,206-04 

6 

-8 , 346-04 

-3.42E-04 

-2.05E-04 

-1 .31E-04 

1,416-04 

5 

-1 .86E-04 

•1,236-Oa 

-7.75E-05 

-3,226-05 

1,516-04 

4 

-4.68E-05 

-5.42E-08 

9,936-06 

3,196-05 

1 , 4 8E- 0 4 

3 

-1.906-05 

1 ,105-05 

2,656-05 

4,706-05 

1 .22E-04 

2 

-8 , 35E-06 

1,166-05 

2.296-05 

3.59E-05 

7,916-05 

1 

-2,38E-0fe 

2 .0SE-06 

6,076-06 

1 .03E-O5 

2,1 3E-05 

1 I 

= 1 

2 

3 

4 

5 

, f\ 

t - 

5.23E-U 

9 , 63E- 1 1 

7,476-1 1 

-2,616-10 

- 1 , 066-09 

9 

4,016—11 

6.96E-! 1 

1,046=10 

=1,406=10 

_ + + rt r _ i n 

- i , Jvi_- » V 

8 

1 ,856-1 1 

2,976-11 

7,876-11 

-5.85E-1 1 

-4.59E-1 0 

7 

-4,61£«12 

-6,436-12 

3,886-11 

-1,326-11 

-2,556-10 

6 

-2.11E-U 

-2.86E-1 1 

2.776-12 

7.83E-13 

-1 .26E-10 

5 

-2.28E-11 

-3.65E-U 

-1 .97E-1 1 

-1.90E-12 

-5,616-11 

4 

-1,656-11 

-3.06E-11 

-2,536-11 

-6,86E-12 

-2.48E-11 


3 "7 ,356=12 
2 -6,546-13 
1 -5,106-13 


j 

10 

9 

8 

7 

6 

5 

4 

3 

2 

1 


1= 1 
8,776-04 
1.08E-OJ 
1.276-03 
1,246-03 
8,346-04 
1,866-04 
4,685-0? 
1,906-05 
8.355-06 
2.386-06 


>1 ,6b£-l 1 
• 1 , 216-11 
• 1 . 37 E-U 


2 . 456-04 
2 , 746-04 
3 . 23 E -04 
3 , 996-04 
S. 92 E -04 
1 , 23 E -04 
5 . 42 E -06 
■l.tOE-OS 
•1 , 16 E- 0 B 
• 2 , 056-06 


•1.686-11 
-1,486-11 
-2,706-1 1 


7 , 186-04 
5 , 946-04 
4 , 486-04 
3 . 21 E -04 
2 , 05 E« 0'4 
7 , 75 E- 0 S 
- 9 . 93 E -06 
- 2 , 656-05 
- 2 . 29 E -05 
- 6 . 97 E -06 


> 7 , 366*12 
■ 1 . 35 E -11 
• 3 . 52 E -11 


6 . 62 E -04 
5 . 27 E -04 
3 . 84 E-Q 4 
2 . 49 E -04 
' 1 , 31 6-04' 
3 , 22 E-Q 5 
- 3 . 19 E -05 
- 4 , 706=05 
- 3 . 59 E -05 
- 1 . 03 E -05 


■ 1 , 44 E-i 1 
- 2 . 15 E -11 
- 3 , 746-11 


1 . 77 E -04 
1 . 02 E -06 
- 8 , 0 0 E -05 
• 1 . 20 E -04 
-T, 4 16-0 4 
- 1 . 51 E -04 
- 1 , 486-04 
- 1 , 226-04 
- 7 , 916-05 
- 2 , 136-05 


V?i QCITY 

*»* K 

6 

0 , 00 E +00 
O.OGE+OO 
0 , 006+00 
0 , GOE + O 0 
0 , OQE+OO 
O.OOE+OO 
0 , 006+00 
O.OOE+OO 

o.ooe+oo 

O.OOE+OO 

*** K 
6 

- 3 . 76 E -05 
3 . 43 E -04 
5 . 38 E -04 
6 . 32 E -04 
6 , 55 E -04 
8 , 3 IE -04 
5 . 98 E -04 
4 , 926-04 
3 . 07 E -04 
8 , 046-05 
*** K 
6 

- 1 , 96 E — 09 
- 1 . 49 E -09 
- 1 . 04 E -09 
- 6 , 72 E -1 0 
- 4 . 04 E -16 
- 2 . 27 E -10 
- 1 . 40 E -10 
- 1 , 016-10 
- 8 . 48 E -11 
- 8 , 62 E - 1 1 
*** K 
6 

3 , 76 E -05 
- 3 , 43 E -04 
- 5 . 38 E -04 
- 6 . 32 E -04 
- 57555-04 
- 6 . 31 E -04 
- 5 . 98 E -04 
- 4 . 92 E -04 
- 3 , 076-04 
- 8 . 04 E -05 


9 a 

3 1 *** 




3 2 *** 

7 


= 3 **» 

7 


s 4 *** 
7 


**** CHECK2 PHJNT OF WH8 


J 

10 

9 ' 

8 

7 

6 

5 

4 

3 

2 


1= 1 

0 , 00006+00 
Q.OOOOE+60 
o , oo ooE + a o 
0 , 00006+00 
0 , OOOQE+OO 
0 , OOOOE+OO 
O.OQOOE+OO 
0, OOOOE+OO 
o , o 0 fi n 6 + o o 


, OOOOE+OO 
OOOOE+OO 
,00006+00 
.OOOOE+QO 
, OOOOE+OQ 
, OOOOE+OO 
, OOOOE+OO 
, OOOOE+OO 

,'VirnCid'i 


11b 1 XUslO 
3 4 

0 , OOOOE + OO 
0 . OOOOE+OO 
0 , OOOOE + OO 
0 , OOOOE+OO 
0 , OOOOE + OO 
0 , OOOOE+OQ 
0 , OOOOE+OO 
0 , OOOOE + OO 
r> n .\c I r n 


Jl* 1 JL=10 
5 

0 , OOOOE + OO 
0 , OOOOE + OO 
0 , OOOOE + OO 
0,00006+00 
0 , OOOQE + OO 
0 , OOOOE+OO 
0 , OOQOE+OO 
0 , OOOOE + OO 


6 

9 

10 

O.OOE+OO 

0 , 006 + 00 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

0, OQE+OO 

0,006+00 

0,006+00 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

0, OOE + OO 

0,006+00 

O.OOE+OO 

O.OOE+OO 

O.OOE+OO 

0,006+00 

O.OOE+OO 

0 , OOE + QO 

O.OOE+OO 

O.OOE+OO 

0,006+00 

O.OOE+OO 

0,006+00 

0,006+00 

O.OOE+OO 

8 

9 

w 

-4, 00E-04 

-5 , 75E-04 

=5.586=04 

=5 , 55E=05 

=2,53E=04 

-2.66E-04 

1,486-04 

-4.81E-05 

-8.58E-05 

2.62E-04 

7.90E-09 

2,246-05 

3 , 03E-04 

1.50E-04 

8.63E-Q5 

2,946-04 

1 .85E-04 

1,296=04 

2.75E-04 

2, 13E-04 

f ,"696 = 04 

2 ,66£«0 4 

2 , 45E-04 

2.10E-04 

2,326-04 

2,336-04 

2,046-04 

7,766-05 

7,426-05 

5,956-05 

8 

9 

10 

-1 ,206-06 

-2,506-06 

-4,0 i E-08 

=1,186—08 

-2,476-08 

*3, 826=08 

-1,096-08 

-2.31E-08 

-3,556-06 

-9,726-09 

-2,146-08 

-3.35E-08 

-8.43E-09 

-1.97E-08 

-3.22E-08 

-6 , 996-09 

-1.74E-08 

=3, 03E-08 

-5,306=09 

=1,346=08 

=2,47E»08 

-3,4 IE-09 

=6, 55E-09 

-1.07E-08 

-1 .66E-09 

4 , 20E=09 

1.62E-08 

-5,476-10 

1 ,42E = 08 

4.17E-08 

6 

9 

10 

3.99E-04 

5, 75E-04 

5 , 58E-04 

5.55E-05 

2,536=04 

2,666=04 

-1,486-04 

4,80E=05 

8,586-05 

•2 , 62E-04 

-7.90E-05 

-2.25E-05 

«'3,‘<J3r=04 

-T750r«04 

-8.64E-05 

-2.94E-04 

•1,85E=04 

•1.29E-04 

-2 , 75E-04 

=2, 136 = 04 

-1,696=04 

-2, 66E-04 

•2,456=04 

•2.10E-04 

-2.32E-04 

-2.33E-04 

•2,046=04 

•7,766=05 

-7.42E-05 

-5,946=05 

. • • ■ f 

a 

' 9 ~ 

10 


Table 3, 
Anqular 


o, OOOOE+OO 
0 , OOOOE + OO 
0 , OOOOE + OO 
0 , OOOOE + OO 
0 , OOOOE + OO 
0 , OOOOE + OO 
0 , 0 QOOE + OU 
0 , OOOOE + OO 


U.WVVVtTVU 

0 , OOOOE+OO 
0 , OOOOE + OO 
0 , OOOOE+OO 
0 , OOOOE + OO 
0 , OOOOE + OO 
0 , OOOOE + OO 


0 , OOOOE + OO 
0 , OOOOE+OO 
0,00006+00 
0 , OOOOE+OO 
0 , OOOOE+OO 
0 . OOOOE+OO 


V.UUUU&-UW 

0 , OOOOE + OO 
0 , OOOOE+OO 
0 , OOOOE+OO 
0 ,OOOOE+O 0 
0 , OOOOE + OO 


0, OOOOE + OO 0, 
0, OOOOE+OO 0, 
0,00006+00 0, 
0, OOOOE+OO 0, 
0, OOOOE+OO 0, 
0,00006+00 0. 


8 

Momentum 


,00006+00 
OOOOE+OO 
, OOOOE + OO 
>OOOOE+OQ 
OOOOE+OO 
OOOOE+OO 
OOOOE+OO 
OOOOE+OO 
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Table 3.9 Kinetic Energy of Turbulence 

• k.E, of TU'HSULENCE 








*** K 

3 1 OS* 




J 13 1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

10 

5.92E+00 

1.32E+01 

1.94E+Q1 

1.61E+01 

6.99E+00 

2,586+01 

2 , 33E + 0 1 

2.19E+01 

2.UE + 01 

2.09E+01 

9 

1 ,6i£+0 1 

7 , 9GE+G 1 

9,626+01 

8,396+01 

6 , 75E+0 1 

1 , 03E+02 

8.09E+01 

6,506+01 

5 , 42E + 0 1 

4.71E+01 

8 

2.076+01 

1,136+02 

1.29E+02 

1.12E+02 

8,836+01 

1,196+02 

9,936+01 

8.13E+Q1 

6,816+01 

5,896+01 

7 

3,235+01 

1,346+02 

1.48E+02 

1,306+02 

9,956+01 

1,216+02 

1.07E+02 

8 , 95E+0 1 

7.58E+01 

6.60E+01 

6 

6,195+01 

1.38E+02 

1.49E+02 

1,346+02 

1.02E+02 

1, 15E+02 

1. 076 + 02 

9,196+01 

7.91E+01 

6, 9bE+0 1 

5 

5,486+01 

9.20E+01 

1 , 1 5E+02 

1,156+02 

9.24E+01 

1,056+02 

1.016+02 

8 , 90E + 0 1 

7.83E+01 

7 , 00E+0 1 

4 

1,446+01 

3.24E+01 

5,905+01 

7.64E+01 

7.32E+01 

9.22E+01 

8.96E+01 

8.17E+01 

7,416+01 

6.776+01 

3 

2 , 0 7E + 0 1 

2,496+01 

3.39E+01 

4,686+01 

5,436+01 

7,706+01 

7 , 60E + 0 1 

7.23E+01 

6,806+01 

6.38E+01 

2 

4,015+01 

4,295+01 

4.15E+01 

4,206+01 

4.60E+01 

6.53E+01 

6.57E+01 

6, 446+01 

6.24E+01 

6 ,006+0 1 

1 

2, 1 7E+0 l 

4.08E+01 

4,286+01 

4.10E+01 

4.30E+01 

6.19E+01 

6.12E+01 

6,036+01 

5.92E+01 

5.78E+01 







*** K 

a 2 *** 




J I 

- 1 

2 

3 

4 

5 

6 

7 

8 

9 

to 

10 

5,925+00 

1.32E+01 

1 .94E+01 

1 .61E+01 

6.98E+00 

2.58E+01 

2,336+01 

2.19E+01 

2 , 1 1 E+0 1 

2.10E+01 

9 

1,825+01 

7 , 9QE+ 0 1 

9,626+01 

8,396+01 

6,756+01 

1.03E+O2 

8,106+01 

6.50E+01 

5.42E+01 

4 , 7 1 E+0 1 

8 

2, 0 7E + 0 1 

1,135+02 

1 .29E+02 

1.12E+02 

8.83E+01 

t , 19E+02 

9.93E+01 

8,136+01 

6 ,8 1E + 0 1 

5,896+01 

7 

3,235+01 

1.34E+02 

1.48E+02 

1.30E+02 

9,956+01 

1.21E+02 

1,076+02 

8,956+01 

7.58E+01 

6,606+01 

8 

8,195+01 

t,3SE+02 

1 , 496+02 

1 s 34E+02 

1 .02E+02 

1.15E+02 

1.076+02 

9.19E+01 

7,916+01 

6 , 96E+0 1 

e 

5.48E+0 1 

9 ,20E + 0 1 

1 , 15E+02 

1,156+02 

9.24E+Q1 

1,056+02 

1,016+02 

8.90E + 0 1 

7,836+01 

7 , OOE + O i 

4 

1,445+01 

3.24E+01 

5.90E+01 

7,646+01 

7,326+01 

9.22E+01 

8.96E+01 

8 , i 76+0 i 

7,416+01 

6,776+01 

3 

2,07E+01 

2.49E+01 

3.395+01 

4.68E+01 

5.43E+01 

7.70E+01 

7,606+01 

7.23E+01 

6,806+01 

6.38E+01 

2 

4,015+01 

4.29E+01 

4, 155+01 

4.20E+01 

4.60E+01 

6.53E+01 

6.57E+01 

6.44E+01 

6.24E+01 

6 , OOE + O t 

1 

2,175+01 

4,086+01 

4 , 26E + Q 1 

U.10E+01 

4,306+01 

6,1 96+0 1 

6.12E+01 

6,036+01 

5 , 92E + 0 1 

5. 786+01 







*** K 

a 3 *** 




J I 

= 1 

2 

l 

4 

5 

6 

7 

8 

9 

10 

10 

5.92E+00 

1 ,32E+01 

1 .94E+01 

1.61E+01 

6 , 98E+00 

2.58E+01 

2.33E+01 

2,196+01 

2.11E+01 

2,106+01 

9 

1,826+01 

7 , 90E+0 1 

9.62E+01 

8.39E+01 

6.75E+01 

1 .03E+02 

8,106+01 

6.50E+01 

5.42E+01 

4,716+01 

8 

2,075+01 

1,136+02 

1 ,29E+02 

1.12E+02 

8.83E+01 

1,196+02 

9.93E+01 

8.13E+01 

6.81E+01 

5,896+01 

7 

3,235+01 

1,346+02 

1 , 46E+02 

1,306+02 

9,956+01 

1.21E+02 

1 , 07E + 02 

8.95E+01 

7, 58E+0 1 

6,606+01 

6 

6, 19E + 0 1 

1.38e+02 

1 .49E+02 

1 .34E+02 

1.02E+02 

1.15E+02 

1 , 07E + 02 

9.19E+01 

7.91E+01 

6,966+01 

5 

5,485+01 

9 , 20E + 0 1 

1.13E+02 

1.1SE+02 

9.246+01 

1.05E+02 

1 .01E + 02 

8 , 90E + 0 1 

7.83E+01 

7,006+01 

4 

1,446+01 

3.24E+01 

5.90E+01 

7.64E+01 

7.32E+01 

9,226+01 

8 g 96E + 0 1 

8.17E+01 

7.41E+01 

6.776+01 

3 

2.07E+01 

2.49E+01 

3 , 39E + 0 1 

4 , 68E+0 1 

5.43E+01 

7.70E+01 

7.60E+01 

7.23E+01 

6.80E+01 

6.38E+01 

2 

4, 0 1E + 0 i 

4,296+01 

4.15E+Q1 

4 , 20E+0 1 

4,606+01 

6,536+01 

6.57E+01 

6.44E+01 

6,246+01 

6, OOE + O 1 

1 

2,175+01 

4 , 08E+0 1 

4.26E+01 

4.1QE+01 

4,306+01 

6,196+01 

6,126+01 

6.03E+01 

5,926+01 

5,786+01 







*** K 

a 4 *** 




If 

M 

T 

3 

3 

4 

5 

6 

7 

8 

9 

10 

10 

5,926+00 

1,326+01 

1.94E+01 

1.61E+01 

6.99E+Q0 

2,586+01 

2.33E+01 

2.19E+01 

2.11E+01 

2.09E+01 

9 

1 , 82E+0 l 

7.90E+01 

9,626+01 

8 , 39E+0 1 

6 , 75E+0 1 

1 ,036+02 

8,096+01 

6.50E+01 

5.42E+01 

4.71E+01 

8 

2,076+01 

1,116+08 

1.29E+02 

1.12E+02 

8.83E+01 

1,196+02 

9.93E+01 

8, 13E+01 

6,8 1E+0 1 

5.89E+01 

7 

3, 235+01 

1, 346+02 

1.48E+02 

1.30E+02 

9.95E+01 

1,216+02 

1 , 07E+02 

8,9SE+01 

7.58E+01 

6,606+01 

6 

6, 19E+0 1 

l,3«E+08 

1,496+02 

1 , 34E+02 

1,026+02 

1 , 15E+02 

1.07E+02 

9.19E+01 

7.91E+01 

6.96E+01 

5 

5,48E+0i 

9,206+01 

1,156+02 

1,156+02 

9, 24E+0 1 

1,056+02 

1.01E+02 

8.90E+01 

7.83E+01 

7.006+01 

4 

1,446+01 

3.24E+0! 

5.90E+01 

7.64E+0 l 

7.32E+01 

9,226+01 

6,96E+0r 

87 17 E+ 01 

77ne>or 

6,776+01 

3 

2, Of E+0 1 

2.496+01 

3,396+01 

4,686+01 

5.43E+01 

7.70E+01 

7.606+01 

7 ,23E+01 

6.80E+0 1 

6.386+01 

2 

4,016 + 0 1 

4,296+01 

4,i5E+0i 

4,206+01 


6,536+01 

6,576+01 

6,44E+01 

6,246+01 

6.Q0E+01 

1 

2, 17E+0 i 

4.08E+01 

4,266+01 

4.10E+01 


6,196+01 

6, 12E+01 

6, 03E+0 1 

5, 92E+0 l 

5.76E+01 
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i a n l i i 

I uD 1 6 O . iU nuSO i Li L0 


i viiiifJC T G v'J ! C 


• TEMPERATURE • * 

*** K a 1 *** 


J I 

= 1 

a 

3 

4 

5 

6 

7 

S 

9 


10 

10 

5,566+0? 

3.526+02 

5.52E+02 

5,686+02 

6,396+02 

1,396+03 

1,446+03 

1.46E+03 

1.48E+03 

1 

,506 + 03 

9 

5,516+02 

5.41E+02 

5.40E+02 

5,546+02 

6,336+02 

1.41E+03 

1.50E+03 

1.53E+03 

1.56E+03 

1 

.5BE + 03 

8 

5,476+02 

5,356+02 

5.34E+02 

5,466+02 

6,226+02 

1 , 46E+03 

1.56E+03 

1 , 60E+03 

1.62E+03 

1 

, 65E + 03 

7 

5.40E+02 

5,306+02 

5,296+02 

5.39E+02 

6,096+02 

1 .54E+03 

1.65E+03 

1.68E+03 

1.71E+03 

i 

, 73E + 03 

6 

5,276+02 

5,236+02 

5.246+02 

5,336+02 

5 , 97E+02 

1 .66E+03 

1.76E+03 

1.79E+03 

1.80E+03 

1 

.82E+03 

5 

5,076+0? 

5.136+02 

5,186+02 

5 , 29E+02 

5,876+02 

1.79E+03 

1 , 88E + 03 

1.90E+03 

1.91E+03 

1 

, 9 IE + 03 

a 

5.02E+02 

5,086+02 

5, 18E+02 

5,306+02 

5 , 82E+02 

1 .90E+03 

1,956+03 

1 .95E+03 

1.95E+03 

1 

.966+03 

3 

5 , 1 5E+02 

5,286+02 

5.376+02 

5,456+02 

5,846+02 

1.87E+03 

1,896+03 

1.90E+03 

1 ,926+03 

1 

.93E+Q3 

2 

6,076+02 

5 , 98E+02 

5.89E+02 

5,796+02 

5.96E+02 

1,766+03 

1,816+03 

1.85E+03 

1,886+03 

1 

,906 + 03 
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Table 3. 12 Convective Fluxes Through the Combustion Chamber Crossections 



SECTION 4 

PARAMETRIC STUDIES OF INJECTOR ANOMALIES 


4.1 OBJECTIVE 

The objective of the present study is to demonstrate a capability to simulate 
two-phase spray flow, heat transfer, evaporation and combustion within the 
rocket engines for predicting the effects of reactant injections nonuniformity 
upon local and overall heat transfer. For this purpose, six test cases have 
been considered. The basic test case has been set up based on the data 
specified by NASA. First, four model sensitivity tests were performed to 
examine the influence of physical and numerical parameters on the fluid flow 
and combustion. The final sixth test case considered an injection nonuniformity, 
in which 25% of the central part of the fuel injector had been blocked. 

All results are presented in graphical form. 

4 . 2 THE COMBUSTION CHAMBER AND INJECTOR GEOMETRY 

The configuration of the model combustion chamber employed in the present 

calculations is shown in Figure 4.1. The combustion chamber is 0.33m long. 

Chamber diameter at the injection plane is d = 121.92mm and at the throat 

o 

d^ = 66.04mm, Contraction ratio is approximately 4. 

Figure 4.2 presents details of the fuel and oxidizer injector plane. Fuel is 
supplied to the injector through the fuel connector (marked FUL on figure 4.2) 
and axially enters the triplet injectors. Liquid oxygen enters radially 
(connector OX on Figure 4.2) and is distributed through the annular collector 
to all oxygen triplet inlets. 

Detail C on Figure 4.2 presents the L0X-RP1-L0X triplet. The triplet 
arrangement within the injector plane is also shown in Figure 4.2. Note that 
due to the specific hole arrangement the injection plane has two symmetry lines. 
Therefore, a 90° sector has been used for the numerical simulation studies. 

The grid arrangement for the 90° sector of the combustor is shown in Figure 4.3. 

A uniform nonorthogonal grid in x-axial, y-radial and z(=0) circumferential 
directions with NX : N Y : NZ - (21:8:4) is employed for all test cases. Figure 4.3 
presents the injector hole distribution and the y-0 grid. 

The calculation domain in the axial direction extends up to the throat section 
of the nozzle. 
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4 • 3 BASIC TEST CASE SPECIFICATION 

A fully coupled two-phase spray flow, evaporation heat transfer and combustion 
with a uniform reactant injection distribution has been considered. Table 

4.1 presents a summary of input data specifications for the basic test case. 

It includes: 

- geometrical data, 

- physical parameters, and 

- injector data. 

During the calculations it has been assumed that oxygen enters the combustion 
chamber in a fully vaporized form, while kerosene (RPl) enters in the form 
of discrete jets of the liquid fuel sprays. To describe the gas phase motion 
(oxygen + combustion products) Eulerian frame of coordinates has been used. 

The liquid phase (fuel) is described in Lagrangian coordinates. The spray 
is represented by 80 individual droplet parcels (16 slot and 5 droplet 
diameters). The liquid spray jet, at each injection slot, is represented by 5 
droplet sizes D . 

Figure 4.4 presents assumed droplet distribution function and selected 
droplet diameters (D m = 5, 10, 20, 30, 40 and 60ym). 



Figure 4.4 Droplet Distribution Function 
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Table 4.1 Basic Test Case Data Specification for 
Combustion in Rocket Engine Calculations 


SPECIFICATION 

DATA 

UNIT 

Chamber radius at injection plane 

60.36 

mm 

Chamber radius at the throat 

33.02 

mm 

Chamber length (to the throat) 

330 

mm 

Injector type UNLIKE TRIPLET 

284 

— 

Total number of triplets 

37 

-- 

Number of fuel holes in 90° sector 

9.25 

— 

Number of oxygen holes in 90° sector 

19.0 

— 

Number of grids (NX:NY:NZ) 

21:8:4 

— 

Total fuel flow rate 
through 90° sector 

4.632 

1 bm/sec 

(0.525) 

kg/sec 

Total oxygen flow rate 

13.924 

1 bm/sec 

through 90° sector 

(1.578) 

kg/sec 

Fuel type ( c n,g6 H 23.23^ 

RP1 

— 

Fuel molecular weight 

167 

g/mol 

Fuel boiling temperature (RP1) 

487 

°K 

Heat of Vaporization (RP1) 

56.5 

kcal / kg 

Heat of Combustion (RP1) 

236700 

(10457) 

J/kg 

( kcal / kg) 

Liquid Fuel Inlet Temperature 

60 

(288.6) 

OF 

(°K) 

Density of Liquid Fuel at 300 °K 

800 

kg/m 3 

Laminar Vapor Diffusion Coefficient 

1.5 • 10" 4 

kg /ms 

Vapor Specific Heat 

2000 

J/kg K 

Oxygen Inlet Temperature 

150(83) 

R (K) 

(assumed vapor oxygen temperature) 

150 

K 

Operating Pressure 

41.37 • 10 5 

N/m 2 

(600) 

(psia) 

Droplet Median Diameter 

27 

ym 

(for distribution function see 
Figure 4.5) 









The mass-median droplet diameter of the droplet distribution function 
(Dpm s 20-30pm) for the triplet injector has been selected based on experimental 
data of FANG (Reference 1.). For jet diameter 1.85 mm the approximate Dpm is 
27um (See Figure 4.5). All droplets are assumed to be injected axailly from the 
liquid fuel nozzle. 

4 • 4 TEST CASES FOR THE PARAMETRIC EVALUATION. 

The intent of the model sensitivity study was to asses the behavior of the 
model when flow conditions, empirical factors and numerical parameters were 
altered about selected mean values: Under the current contract this task was 

limited to five test cases. The selected five test cases included variation in: 

- droplet evaporation rate formula, 

- vapor fuel reaction rate formula, 

- droplet diameter distribution function, and 

- liquid fuel injection distribution. 

Table 4.2 presents the basic test case (TO) default parameters and five 
test cases (T1 to T5) parameter specification summary. Note that in test 
cases T1 - T5 only one parameter at a time has been altered; the remaining 

parameters were specified as in the basic test TO. In all cases a 21: 8: 4 

grid has been used. 

y Table 4.2 Test Cases for Parametric Evaluation 


TEST 

NUMBER 

TEST CASE DESCRIPTION 

DIFF. COEF 

(M) 

'ms' 

REACTION RATE 

CONST. P CM 
rU 

l OF FUEL HOLES 
BLOCKAGE 

TO 

Base case 

1.5 • 10" 4 

io 10 

0 

T1 

Doubled evaporation rate 

3.0 • 10" 4 

io 10 

0 

T2 

Halved evaporation rate 

1 

o 

r — 1 

• 

LO 

• 

o 

10 10 

0 

T3 

Increased reaction rate 

1.5 • 10" 4 

10 12 

0 

T4 

Altered droplet distri- 
bution 

1.5 • 10" 4 

io 10 

0 

T5 

Blocked fuel holes 
near axis 

1.5 • 10~ 4 

io 10 

~25 


Further details of the test cases are described below: 
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Mass-Median Drop Radius 



m 









Influence of Diffusion Coefficient in the Evaporation Rate Formula 

First two test cases T1 and T2 were designed to investigate the influence 
of the vapor diffusion coefficient of V on the rate of evaporation of the 
liquid droplets. The evaporation rate formula is exressed as: 


dD = 
dt 


4pP . 
C D 1 
P 


(1+B) 


where 


C 4 (T-Tsat) 
L 


(4.1) 


and in the basic test case, the diffusion coefficients V has been estimated 
“62 

to be: 10" m /s (see reference 2 page 204). The pV is obtained as: 


n f) Hi H . t) ~ 41 . x 10 • 196 ~ (1 ~ ? 1 x ia"^ kg (a?) 

pV 8314-T..X V ~ 8314*600 ° 1U ~ ”** x 10 me (4.2) 


sat 


ms 


In the basic test case (TO) value of pV equal to 1.5*10 ^ —2. has been used. 


For model sensitivity studies, this value has been respectively doubled (Test 
Case Tl) and reduced to half (Test Case T2). 


Influence of the Reaction Rate Constant 


The two step combustion reaction rate scheme employed in the model assumes 
Arrhenius reaction rate expressions in the following form: 

R fu ■ R fu<"fu> afU ('"ox> • exp(-E fu /RT) 


R co = R co< m co> aC ° <V> bc °<% 2 > CC ° ex P(- E co /RT > 


(4.3) 

(4.4) 


where a f(j , b fu# a co , b co , c co = 1 and 
in i 

P fu = 1.10 iU (kg/m s) E fu = 18000 (cal/g) 

P = 7.10 10 E = 1860 

CO CO 

The fuel reaction constant P^ u in the reaction mechanism is less reliable. At the 
same time this reaction provides most of the hear release. In the present 
study a test case with P^ u increased to 1.10 (two order increment) has 
been considered to study changes in the flame structur. 
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Influence of the Droplet Distribution Function 

Liquid fuel injected into the combustion chamber through the axial hole of the 
"triplet L0X-RP1-L0X" injector is atomized and sprayed in the form of droplets 
of different diameters. Exact specification of the spray characteristics 
for the chamber operating conditions is very difficult. Approximate 
specification of the droplet diameter distribution function at the injection 
location requires knowledge of median droplet diameter (D m )» and minimum and 
maximun diameters ( D m n - n , D max ). Experimental data by Fang [1] indicate that 
the median droplet diameter for the triplet injector and RP1 fuel can be taken 
equal to 27 ym„ The maximum droplet diameter is estimated as 80 -120 pm. 

Figure 4.6 presents two assumed droplet distribution functions 
at the injection point for basic test case TP (solid line) and for test 
case T4 (dotted line). In both cases median droplet diameter is similar. In 
the test case T4 maximum droplet diameter has been reduced and number 
of smaller diameter droplets (15 to 30ym) increased. 

Table 4.3 presents discrete droplet distribution function (D - droplet 
diameter and P-droplet population number) for both test cases. 



Figure 4.6 Droplet Distribution Functions at the Inlet (Injection) 

to the Chamber for Test Cases TO - T4. 
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Table 4.3 Droplet Distribution Function For Basic Test Case TO and Test Case T4. 


- — 
TO - TEST CASE 

T4 - TEST CASE ; 

D 

P 

D 

P 

0 


0 


pm 

.... 

ym 

- 

5 

20 

5 

10 

10 

35 

10 

30 

20 

25 

20 

40 

40 

15 

35 

15 

60 

5 

50 

5 


Influence of Injection Nonuniformity 

The injector nonuniformity, in general, would require specification of the 
fol lowing: 

a) location of partially or fully blocked fuel holes, 

b) location of partially or fully blocked oxygen holes, and 

c) percentages of blockages. 

During the present study, only one test case (T5) of the injection nonuni formity has 
been considered for the computations. 

It has been assumed that 25% of the fuel holes have been blocked. The blocked 
holes are located symmetrical ly near the chamber axis so the calculation 
domain stays unchanged as in the basic case TO. 

Figure 4.7 presents the injector hole arrangement, fuel holes blockage 
and (y-0) grid distribution at the injector plane. The blocked fuel entry area 
comprises 24.3% of the nominal fuel entry area. The fuel mass flow rate through open 
holes has increased so as the total flow rate is maintained to be the same 
as in basic test case TO. 

4,5 RESULTS OF BASE TEST CASE 
Presentation of the Results 

Figure 4.8 presents typical velocity and pressure contours on an axial (XY) ~ 
plane of the combustor. Figures 4.9 and 4.10 show contours of absolute temperature 
(in °K) and contours of vaporized fuel concentrations in all four (X-Y) planes 
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injector plane 


Blocked 
Fuel Holes 


K»4 



Figure 4.7 Injector Plane with Blocked Fuel Holes 
Near the Axis of the Chamber (Test Case T6). 

Total No. of Holes in 90° Sector: 

Blocked No. of Holes 
% Blockage 


9i 

2 X 

24.3% 
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XV PLANE 1 
VELOCITY PLOTS 
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9 — 1 - 116E+05 
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Figure 4.8 
Plane (XY) 


Velocity Vectors and Pressure Contours on an Axi 
of the Combustor. (Basic Test Case Tl). 
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Figure 4.9 Temperature Contors in Axial (XY) Planes; Bas 
Test Case (Tl) . 
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Figure 4.10 Contours of Gaseous Fuel Mass Fraction in XY Planes. 
(Basic Case Tl). 


of the combustor. Figure 4,11 presents fuel concentrations at three 
selected cross-sections of the combustor. Figure 4.12 shows the distribution 
of the convective heat transfer coefficient on one (X-Z) plane along the 
combustor wall (developed cylindrical surface). Figure 4.13 presents the 
heat transfer coefficients variation along the cylindrical wall for each 
0 - plane. Some of the results of the basic test case TO, such as total 
gaseous fuel flow rate, liquid spray rate, etc, are presented and discussed in 
the section describing parametric evaluation. 

Discussion of the Results 

The flow field represented by the velocity vecotor and pressure contours 
is typical of nozzle flows. Largest pressure gradient exists near the 
combustion throat, and the lowest pressure is also located in the same region. 

In the front of the combustor, the gas velocities are relatively low. At 
a distance approximately equal to the chamber radius, a large increase in 
velocities takes place. At that distance the flame fronts of the individual 
triplet injectors are joined and large temperature gradients occur (see figure 
4-9). 

Figure 4.9 indicates that the flame front is nonuniform in both circumferential 
and axial directions. Temperature rise takes place in shorter axial distance 
at mid-radii than at small radii (i.e. near the axis). Note that at 2nd and 
3rd XY - plane near the chamber wall there is a recirculation region carring 
hot combustion products towards the injection plane. 

The largest fuel concentration gradients exist in the front part of the chamber 
where droplets enter hot combustion zone and intensive evaporation takes place. 

The smallest droplets evaporate in a short distance after entering the reaction 

zone creating a region of significant interphase mass transfer. The largest 

droplets penetrate further into the chamber and create a zone of reacting 

vapor fuel along the trajectory. Figure 4.14 presents qualitative comparison 

between the experimental results (photograph of single triplet flame [1] and calculated 

gaseous fuel concentration profiles at K = 2 plane. In both diagrams (Figure 

4-14), an enlongated flame shape can be seen. Note that the photograph 

represents only a "hot luminous" zone where unburned hydrocarbon creates a 

luminous zone. No conclusion can be drawn regarding high temperature zone from 

this photograph. Figures 4.10 and 4.11 present the gaseous fuel mass fraction 
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jgure 4.12 Heat Transfer Coefficient on the Cylindrical Wall 
(Development of X-Z Plane). Case T1 
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Figure 4.13 Distribution of convedtive heat fluxes along the cylindrical wall 



HIGH-SPEED PHOTOGRAPH OF OFO TRIPLET 
CARBON FORMATION PHENOMENA, TEST 116 



4.14 Comparison between predicted fuel mass function contours and 
experital flame photograph. 
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contours at different cross-sections of the combustor. Long wakes of the 
unburned vapor fuel are generated along the liquid jet trajectories (see figure 
4.10). In Figure 4.11, the streams of unburned vapor fuel are represented 
by concentric semicircles. 

Distribution of heat transfer coefficients r at the cylindrical wall is 
shown in Figure 4.12 (contour map) and in Figure 4.13 (distribution curves). 

From Figure 4.13, it is seen that the heat transfer coefficient is increasing 
with the distance, with the maximum being near the throat region. Similar 
observation has been reported in experiments. Note that r is quite nonuniform 
in the circumferential direction, especially in the upstream part of the 
chamber. The heat transfer curves at k = 2, 3, and 4 are quite similar. At 
k = 1, however, heat transfer coefficients are significantly smaller (see also 
Figure 4.12). 

The outline of the grid and injection hole arrangement, shown in Figure 4 J. 2 , 
reveals that as compared with k = 2., 3, and 4 in the k = 1 plane, oxygen 
holes are at the shortest distance from the wall. This results in lower 
velocity and temperature (cooling effect) and lower heat transfer coefficient H 

The cooling effect may have been overpredicted due to grid coarseness at 
the large radii. Additionally, the influence of the oxygen jet is "averaged" 

(or smothered) over the entire grid volume. For better resolution, calculations 
with finer grid or alternative injection specification would be required. 

4 . 6 RESULTS OF PARAMETRIC STUDIE S 

TEST T1 and T2: Influence of Diffusion Coefficient V in Evaporation Formula 

The evaporation rate dD/dt is proportional to the diffusion coefficient V 

(see Equation 4.1). Assumed variants in the diffusion coefficient, viz: 

T1 : pV = 3. Ox 10~ 4 kg/ms 

TO : pV ~ 1.5x 10 -4 kg/ms 
-4 

T2 : pV = 0.75x 10 kg. ms 

is approximately equal to the limits of uncertainty for this factor. 

Figure 4.15 presents the axial variations of liquid fuel flow rate calculated 
for test cases TO, Tl, and T2. The results show that: 

a) larger V creates higher evaporation rates and shorter penetration 
of liquid fuel jets; 
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b) in case Tl, liquid is totally evaporated within 3/4 of the chamber 
length; 

c) in the other two cases, only droplets with the largest initial 
diameter D q can be found at the throat. The liquid fuel flow 
rates at the exit are: 


for Tl 

: *L1, = °- 

. {i.e. OS of rfi L1q _ 1nJ) . 

for TO 

: = 3.10- 4 

• O.e. 0.01* of iti Llq _ jnj h 

for T2 

! *L1q \ 4 - 10 " 3 

. (i.e. 0.7* of H L1q _ 1nj ). 


Detailed flow patterns indicate no major differences between the flow field 
and heat transfer in all three cases. 

Test T3: Influence of the Fuel Reaction Rate Constant Ppy 

Figure 4.16 presents the axial variations of gaseous fuel flow rate for two 
different fuel reaction rate constants; viz: 

10 1 

Case TO with Ppy = 10 kg/rn s ; and 
Case T3 with Ppy = 10^ kg/tn^s. 

In the case of larger Ppy (Case T3) vapor fuel flow rates in the front part 
of the chamber are singificantly smaller that those in the basic case TO. From 
temperature contours (Figures 4.9 and 4.17) it can be seen that in the case 
of larger Ppy, the flame front is located closer to the injector plane. This 
seems physically plausible. 

Test T4: Infulence of Alternative Droplet Distribution Function 

By changing the droplet distribution, there are no significant changes in the 
calculated flow field and heat transfer characteristics. In Test T4, flame 
front is slightly closer to the injector plane. Inspection of the droplet 
distribution functions (Figure 4.6) reveals that in Test Case T4 there is a 
smaller amount of the finest drops (0-10ym) and therefore less vapor fuel 
will be evaporated in the front part of the combustor. This probably creates 
mixture conditions in Case T4 closer to the stoichiometric than that in base 
case TO. 

Test T5: Influence of Injector Nonuniformity (Blocked Fuel Holes Near Axis) 

Figure 4.7 presents injector plane with four blocked fuel holes in the 90° 
sector. Blocked fuel entry area is 24.3% of the nominal fuel entry area. As 
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X/L 

4.15 Axial Variation of Liquid Fuel Flow Rate (for 45° Sector of the Chamber); 

for (Tests TO, Tl, and T2). 



4.16 Axial variation of Gaseous Fuel Flow Rate (for 45° Sector of the Chamber); 

for (Tests TO and T3) . 
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Figure 4^17 Fuel Mass Fraction Contours Test Case T4-Higher Fuel Reaction Rate Constant 
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Figure 4.18 Temperature Contours in XY-Planes. Test Case T5- Near Axis Blockage of 

25% of the Fuel Entry. 
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Figure 4„ 19 Heat Transfer 



icient on the Cylindical Wall Test Case T6- 
ocked 25% of the Fuel Entry. 



c / rcumfe re n flat 




compared to Cases T1 to T4, this test case (T5) shows much more sirtgificant 
differences in the heat transfer and combustion within the chamber. Figure 
4.18 presents the temperature contours (°K) in four axial plans. The central, 
oxygen rich core is significantly cooler than the annular, fuel rich, zone. 

Steep gradients in the radial direction in the middle of the combustor 
represent the enlongated flame front. 

Comparison of Figure 4.12 and 4.19 indicates that convective heat transfer 
coefficients at the combustor wall are nearly the same for the basic case TO 
and blocked fuel entry case T5. This is due to the "blanketing" of the cylindrical 
wall by the hot reactive streams of operating fuel spray jets. An equivalent 
blockage near the combustor wall, rather than near the axis, is likely 
to have a much larger influence on the "wall" heat transfer coefficients. 

Further parametric studies, were not included in the current work scope, and are 
recommended for future. 
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SECTION 5 

MODEL IMPROVEMENTS AND RESEARCH NEEDS 


The applications reported in Sections 3 and 4 of this report have demonstrated 
the basic capability of simulating three-dimensional, two-phase flows with 
evaporation and combustion. The use of Langrangian technique for liquid fuel 
drops has provided necessary flexibility of predicting effects of small 
changes in injector plate geometry and/or flow conditions. Results of 
parametric studies have shown that the numerical model responds to changes 
in physical parameters in a plausible manner, and thus it can be utilized for 
the identification of sensitive parameters. 

Several improvements and studies can and should be made for the code to be 
a valuable tool in design and performance analysis of rocket engines. The 
improvements are desirable for both mathematical models of physical processes 
built into the code as well as for numerical methods. and solution algorithm 
employed. The code also requires verification and sensitivity study for both 
single and two-phase flows with and without mass transfer and combustion. 

The following three subsections describe recommendations for (a) physical 
model improvements; (b) numerical model improvements and (c) code validation 
and verification studies. 

5 . 1 IMPROVEMENTS IN MATHEMATICAL MODELS OF PHYSICAL PROCESSES 

There are several possible improvements in simulating the spray flow, evaporation 
and combustion processes implemented in current version of the code. Further 
refinements can be made in simulating the injector and spray atomization. This 
section briefly describes basic ideas of selected possible improvements. 

Spray Atomization in Triplet Injectors 

The unlike impinging elements accomplish mixing and atomization by direct 
impingement of fuel and oxidizer jets. Atomization takes place in the immediate 
vicinity of the impingement point. The effect of oxygen jet impingement on 
the central fuel stream scatters the fuel stream in the direction normal to the 
injection plane (Figure 5.1). 
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O OUTER ORIFICE 

OUTER ORIFICE FLUX 

• INNER ORIFICE 
INNER ORIFICE FLUX 



(a) Inner and outer jets (b) Outer stream momentum 

momentum balanced much greater than inner 

stream momentum 

Figure 5.1 Triplet Injector Spray Mass Distribution 

Experimental data exists (see Review in [3] ) for specifying a range of droplet 
beams at the atomization point. Such information should be utilized in future 
studies and incorporated in if necessary in the code for regular use in analysis 
of injector anomalies. 

Turbulent Droplet Diffusion 

In the present version of the code, turbulent droplet diffusion has been 
neglected. This must be improved for future studies. The task of including 
turbulent particle diffusion into the model requires investigations of various 
possible approaches. One diffulty is that the fluid properties are constant 
within a cell, while the particles are tracked through a cell in a series of 
time steps smaller than the cell dimensions. Thus the particles at several time 
step positions within a cell see only one set of mean flow properties. 

There are two basic approximate techniques for simulating the turbulent 
droplet diffusion. 

a) random walk method (Dukowicz [4] ) ; and 

b) diffusive drift method of Jurewiez[5] and Stock [6]. 

In both techniques, an isotropic turbulence is assumed. The second technique 
seems to be more adequate for "fully coupled" calculations of spray combustion 
modeling, and is recommended for inclusion in REFLAN3D-Spray code. 

Accurate Evaluation of Thermodynamic and Transport Properties 

In the present version of the code, most of the liquid fuel properties 

(C v , T p, . , pt c) have been assumed to be constant. More accurate 
v p.nq sat Hiq 

representation would require sepcifi cation of their functional dependence 
on local pressure and temperature. 
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The gaseous species specific heat formula coefficients (5-order polynomial) 
are valid for the "high temperature range" between 500-3000°k. For the low- 
temperature liquid propellants, two sets of specific heat constants are desired. 

5 . 2 IMPROVEMENTS IN NUMERICAL METHOD AND SOLUTION ALGORITHM 
The computational process of the coupled Eulerian-Lagrangian analysis is highly 
nonlinear. In the present calculations, no special relaxation or linerization 
of interphase mass, momentum and energy transfer source terms has been undertaken. 
Calculations with (11 x 8 x 4) and (21 x 8 x 4) grid caused no numerical 

instabilities. However, trial runs with finer grids in the circumferential 

direction (21 x 8 x 8) showed slow convergence and oscillatory behavior. 

These difficulties must be investigated and remedied before the code can be 
regarded suitable for studies of injector anomalies. 

Relaxation Practices 

Investigations of the relaxation practices and linearization methods are desired. 
For the two-phase flow with mass transfer, possible relaxation practices include: 

a) relaxation of the interphase mass transfer rate ; 

b) relaxation of the average liquid property (f>-| ^ in the transfer 

rate expression; 

c) relaxation of velocities and enthalpies; and 

d) combination of a, b, and c. 

Accuracy Improvements 

The accuracy of the three dimensional calculations in the Eulerian frame can 
be improved by refining the grid in the region of interest and/or by using 
higher order differencing methods. 

The usefulness of the first practice is often limited by the computer storage, 
especially for the reactive flow calculations where twelve or more differential 
equations must be solved. 

The second approach is higher order finite differencing and has attracted more 
attention in the recent publications '17,8,9], One of the most promising 
methods of the accuracy improvement has been developed by the authors [10,11]. The 
new method called "Multiflux Conservative Differencing" (MCD), has been 
successfully applied in a 2- dimensional calculation. Results obtained with 
the coarse grids 10 x 10 and 20 x 20 compare very well with results on 
50 x 50 or finer grids obtained with the other (e.g. upwind differencing) 
methods. It is recommended that: (a) MCD be incorporated in REFLAN3D-Spray Code, 
and (b) comparative studies be performed to evaluate the improvements due to MCD. 
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The second part of the computational algorithm (Lagrangian part) is relatively 
new in comparison with the Eulerian part and the computational experience 
is still very limited. Some basic studies in the area of Eulerian-Lagrangian 
approach would greatly enhance the computational fluid dynamics. New solution 
schemes for the homogenous and heterogenous combustion and fluid flow 
calculations could be explored. 

For the rocket engine calculations, the first improvement should probably b$ made 
in the calculations of the interphase mass source. In the present version of 
the code, the source terms are calculated at the grid cells through which 
the particles are passing (see Figure 5.2a). 



a) current 



b) PROPOSED 


Figure 5.2 Lagrangian Interphase Source Term Calculations for 
Eulerian Transport Equations. 

Therefore, the mass release is "lumped" in the source at point "P" with no 
source for point N. In reality, however, the droplet trajectory represents 
an average droplet path and one could expect to have the mass transfer to 
the N cell as well. The practice in Figure 5.2b presents an alternate 
approach for the source term calculations around the average particle path. 

The mass transfer from the liquid to the gaseous phase will be distributed to 
both N and P points based on 

a) distance weighting factors or, 

b) volume weighting factors. 

With this practice, implementation of the droplet turbulent diffusion models 
would also be more realistic and economical. Simple test cases on a 2-dimerisional 
grid are recommended as the first step. 
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Improvements in' the Nonorthogonal Grid System 

The nonorthogonal grid system and associated velocity components in the 
present version of the code are shown in figure 5.3a. The orthogonal 
velocity components, u and v are used on the nonorthogonal grid. Within this 
practice it has been assumed that the Up velocity is driven by the (Pp - p w ) 
pressure gradient. This assumption is valid only if the velocity components 
are aligned with the grid lines, or the degree of departure from this condition 
is small (say ^ 30°). A new practice is proposed in which the grid lines 
and velocity resolutes (nonorthogonal) are used figure (5.3b). Most 
of the derivations for preliminary testing have been done under CHAM's in-house 
development project. Its implementation and numerical test in REFLAN3D-SPRAY 
are recommended to enhance the code capability and accuracy. 




Figure 5.3 A Nonorthogonal Grid Systems 
Modification of the Code Structure 

In the present version of the REFLAN3D Code, COMMON statements are used to 
transfer the variables between subroutines themselves. 

A more economical method of passing data to a subprogram is through the 
argument list in the CALL statement. This practice is more economical (small 
storage allocation) and offers greater flexibility by "dynamic storage 
reservation" option. This practice is also preferred in recently specified 
"NASA Standard Requirements for Digital Computer Programs" 

The modification does not require any changes in the bulk of the coding, and is 
recommended for REFLAN3D-SPRAY code. 



Improvements in Solution Algorithm 

The current solution algorithm employs a modified version of the SIMPLE 
procedure [12] for calculating the hydrodynamic field (u-v-w-p). Figure 5-4 
presents the flow chart of the iterative scheme employed. Note that after 
solving the pressure correction equation, velocities are corrected by using 
approximate velocity-pressure coefficients based on truncated momentum equations. 
Therefore, the u, v, and w do not satisfy the full momentum equations with body 
forces, large shear stress and interphase momentum transfer* An alternative 
scheme which can be useful for flow with high shear stresses and/or with 
intensive interphase momentum transfer is shown in Figure 5.5. The scheme 
differs from the SIMPLE algorithm in that after the correction of p, u, v, 
and w, a second correction is performed based on the Poisson equation for 
pressure, obtained from the full (not truncated) momentum equations. 

Velocities are corrected using these secondary pressure corrections and the 
subiterative loop returns to the "continuity conserving" pressure correction 
equation. It is recommended that in the next stage of study, a few variations 
of the proposed scheme be tested and the most economical one employed for 
retention as a permanent feature of REFLAN3D -SPRAY Code. 

5 • 3 VALIDATION AND VERIFICATION STUDY 

The calculations of the three-dimensional two-phase flow in rocket engine 
combustion chambers undertaken in the present project represent the first 
attempt to model 3-D combustion by the Eulerian-Lagrangian algorithm. 

Additional validation and verification tests of the present model are required 
and recommended. The computational studies should be performed for: 

tests for improvements in physical model and numerical algorithm 
(as described in Sections 5.1 and 5.2); 
wider range of physical parameters; 

prediction of single triplet injector flow and combustion 
(Experimental results of the relevant case are available in 
NASA CR 1169006. The geometry of the combustion chamber and the 
photograph of the triplet jet flame are shown in Figure 5.6); 

Injector nonuniformity study for several hole arrangements and 
injection hole blockages. 
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CALCULATE THE INTERPHASE MASS MOMENTUM 
AND ENERGY TRANSFER SOURCE TERMS 



Figure 5.4 Flow Chart, of Current Solution Algorithm 
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Figure 5.5 Flow Chart of Proposed Double Pressure 


Correction Algorithm. 
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j — SINGLE ELEMENT 
INJECTOR' 



Figure 5.6 Experimental Test Chamber Geometry and 0F0 Triplet 
Injection Spray Pattern (Reference 1). 
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APPENDIX A 

INTERPRETATION OF RESULTS AND THE COMPUTER OUTPUT 

The computational results discussed in this report are presented in the form 
of tables of computer output- and graphical plots of velocity vectors and 
contour maps of the dependent Variables. This appendix describes geometry, 
grid distribution notation and interpretation of graphical and tabulated results 
obtained from the REFLAN3D code. 

A . 1 GEOMETRY AND COMPUTATIONAL GRID 

A sample geometry for elliptic flow calculations is presented in figure A.l. 

In this particular case a symmetry plane is assumed and only the 180° sector 
of the flow region is considered as a calculation domain. 

The fluid enters the cylindrical channel through the front plane which is 
partially blocked (Figure A.l), and leaves the channel through the fully 
opened exit plane. 

The selected grid distribution for this geometry is shown in Figures A. 2a and 
A. 2b respectively for two cross-sections viz: 

a) XY - axial plane - aligned with the main stream direction, and 

b) YZ - cross-section plane - perpendicular to the main stream 

direction. 

The coordinate directions are selected as: 

X - axial direction, 

Y - radial direction, and 

Z - circumferential direction, (in radians for polar coordinates). 

The boundries of the calculation domain are identified as: 

West Boundary - located at X = 0, 

East Boundary - located at X = X„,„ 

max 

South Boundary- located at Y = 0 

North Boundary- located at Y = Y 

Low Boundary - located at Z = 0 

High Boundary - located at Z = Z 

s J max 
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Figure A.l 


Geometry for 


test Flow Calculations. 
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Figure A. 2 also presents staggered grid arrangement and the storage location 
of velocity components. Note that "backward boomerang" is used for storing 
all three velocity components viz: U-axial , V-radial and W-circumferential . 

Seperate arrays are allocated for velocities at the last X-face (UEB), last 
y-face (VNB) and last z-face (WHB). All scalar variables (p,T,k,e,m fu . . .) 
are stored in the middle of the grid cell. 

A. 2 INTERPRETATION OF THE PRINTOUT 

Figure A. 3 presents the typical printout of the axial velocity components, 
and appropriate diagrams for the output interpretation. 

All variables (including velocity) are printed in the form of tables; each 
table for a particular axial plane (Z = Const) numbered K=l, K=2, etc. For 
velocity components additional arrays for: 

east boundary axial velocity UEB, 
north boundary radial velocity VNB, and 

high boundary circumferential velocity WHB are printed too. 


Velocities are calculated and printed at appropriate cell faces of the control 
volume. Scalar quantities are printed at the grid centers. All variables are 


printed in SI 
k 
e 
h 
T 
f 
m 
P 
y 


fu 


-units i.e. 

kinetic energy of turbulence 

dissipation rate 

total enthalpy 

absolute temperature 

mixture fractions 

fuel mass fraction 

2 , 


2.2 

m /s 

2. 3 
m /s 

J/kg 

°K 
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pressure (N/m ) relative to the reference point., 
effective viscosity kg/ms 

INTERPRETATION OF THE GRAPHICAL PLOTS 


Graphical representations of the calculation results include: 
velocity vectors in XY, YZ or XZ plane 

isoline contours of scalar variables in XY, YZ or XZ plane. 

The velocity vectors U, are calculated and plotted at the grid centers based 
on the linear interpolation from the velocity components located at the grid 
cell faces (see Figure A. 4). 
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Figure A. 4: Calculation on Velocity Vectors. 


The vector length is taken as: 



Figure A. 5 presents a typical distribution of velocity vectors for the first XY 
plane with an overlayed grid line system and contours of the scalar quantity 
(temperature) on the second XY plane. 

Additional information for velocity vectors includes printout of minimum (VELMIN) 
and maximum (VELMAX) velocity vectors at that plane. 

The legend for the contours of scalar property <|> includes: 
plane orientation eg: XY Plane 2; 

name of the contoured quanity: TEMP; 

minimum and maximum value of <j> within that plane (PHIMAX and PHIMIN); 
and 

contour numbers and appropriate contour levels. 


A- 6 



A-7 


Figure A;5 Example 


XV PLANE i 
VELOCITY PLOTS 
UELN IN 4 . 056E-0 1 

UELMAX £♦ 260EX01 



XV PLANE 2 
TEMP CONTOURS 
PHIMAX PHiniN' 

4.9S19E+02 3 . 0457E+02 

CONTOUR LEVELS ARE 

1 3.2I90E+02 

2 3.3923E+02 

3 3.S65SE-t>02 

4 3.7389E+02 

5 3.9122E+02 

6 4.0855E+02 

7 4.2587E+02 

8 4.4320E+02 

9 4.6053E+02 

10 4 ♦ 7786E+02 



the Graphical Results Interpretation 
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